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ABSTRACT 

We have detected in high resolution spectra of the quasar Q0347-3819 obtained with the 
UVES spectrograph at the VLT/Kueyen telescope over 80 absorption features in the Lyman and 
Werner H2 bands at the redshift of a damped Lya system at z a b s = 3.025. At the same redshift, 
numerous absorption lines of atoms and low ions (H I, D I, C II, C II*, N I, O I, Si II, P II, 
Ar I, Cr II, Fe II, and Zn II) were identified. The z a b s = 3.025 system spans over ~ 80 km s _1 
and exhibits a multicomponent velocity structure in the metal lines. The main component at 
z abs = 3.024855 shows a total H 2 column density N{R 2 ) = (4.10 ± 0.21) x 10 14 cnr 2 and a 
fractional molecular abundance /h 2 = (1.94 ± 0.10) x 10~ 6 derived from the H2 lines arising 
from J = to 5 rotational levels of the ground electronic-vibrational state. In the second strong 
component at z a b s = 3.024684 we were able to set a stringent upper limit to /h 2 of 3.3 x 10~ 9 
(3<7 c.l.). The population of the J levels can be represented by a single excitation temperature 
of T cx — 825 ± 110 K. The ortho:para-H2 ratio equals 3.1 ± 0.3. This ratio is comparable to 
that calculated from hot H2 formation upon grain surfaces (ortho:para-H2 = 3:1), but is higher 
than the 'freeze-out' quantity (ortho:para-H 2 = 1:4) predicted for cosmological epochs at z < 20, 
implying recent formation of the hydrogen molecules following photo-dissociation of pristine H 2 . 
In the main component the properties of the neutral atoms and low ions can be described with 
the same broadening parameter b ~ 3 km s~ x of the H 2 cloud. The rate at which UV photons are 
absorbed in the Lyman and Werner bands is found to be /3q — 2 x 10~ 9 s _1 . Taking into account 
that similar /3q values are observed in the Galactic diffuse molecular clouds, it shows that the UV 
radiation fields in the DLA and in the MW are very much alike. We also found the H2 formation 
rate of i?nn — 3 x 10~ 16 s" 1 at z a b s = 3.025 which is consistent with observations in the Galactic 
disk diffuse clouds where nn 2 is low. Comparing the relative rotational population ratios of 
H2 with the measured C II*/C II ratio, we can infer that nrr — 6 cm -3 in the main component, 
which gives a size of D ~ 14 pc along the line of sight. The mean value of the dust-to-gas ratio 
estimated for the whole z a b s = 3.025 system from the [Cr/Zn] ratio is equal to k = 0.032±0.005 (in 
units of the mean Galactic ISM value) . For the first time we unambiguously reveal a pronounced 
[a-element /iron-peak] enhancement of [0,Si/Zn] = 0.6 ± 0.1 (6a c.l.) at high redshift. The 
simultaneous analysis of metal and hydrogen lines leads to D/H = (3.75 ± 0.25) x 10~ 5 . This 
value is consistent with standard big bang nucleosynthesis if the baryon-to-photon ratio, rj, lies 
within the range 4.37 x lO" 10 < i] < 5.32 x lO" 10 , implying 0.016 < n h h 2 100 < 0.020. 

Subject headings: cosmology: observations — line: identification — quasars: absorption lines — quasars: 
abundances — quasars: individual (Q0347-3819) 
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I. Introduction 

Molecular hydrogen H 2 , being the most abun- 
dant molecular species in the universe, plays a 
central role in star formation processes both at 
low and high rcdshifts (e.g., Shull & Beckwith 
1982). The H2 molecule is known to be an im- 
portant coolant for gravitational collapse of hydro- 
gen clouds at T ~ 1000 K. Fractional abundances 
of H 2 in interstellar H I regions were calculated 
first by Gould & Salpeter (1963) and by Hollen- 
bach, Werner, & Salpeter (1971). They predicted 
the fractional abundance of molecular hydrogen 
/h 2 ^ 10~ 3 , where /h 2 = 2n H J(n H + 2n H2 ) for a 
standard galactic disk cloud with nn — 10 cm~ 3 , 
gas temperature of about 100 K, and the opti- 
cal depth to the cloud center due to dust grains 
r g ~ 0.1. 

The first detection of H2 in the ISM was per- 
formed by Carruthers (1970) with a rocket-borne 
ultraviolet spectrometer. Later, after the launch 
of the Copernicus satellite, the far-UV absorption 
bands of H 2 have been observed in the Galac- 
tic disk toward numerous relatively bright stars 
with visual extinctions Ay of less than about 1 
mag (e.g. Spitzer & Jenkins 1975; Savage et al. 
1977). More recent observations of H 2 with the In- 
terstellar Medium Absorption Profile Spectrograph 
(IMAPS) (Jenkins & Pcimbert 1997), the Orbit- 
ing and Retrievable Far and Extreme Ultraviolet 
Spectrometer (ORFEUS) (Richter et al. 1998; de 
Boer et al. 1998), and the Far Ultraviolet Spec- 
troscopic Explorer (FUSE) (Moos et al. 2000) 
have shown the ubiquity of H 2 in the disk and 
halo of the Galaxy and in the Magellanic Clouds. 
In the FUSE mini-survey, H 2 was detected in the 
sightlines with neutral hydrogen column densities 
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N(H I) ~ (3 6) x 10 20 cm" 2 , the fractional 
abundances /h 2 ^ 7.4 x 10~ 6 3.4 x 10 -2 and 
the excitation temperatures T cx — 180 -j- 540 K 
(Shull et al. 2000). 

In these observations the range of N(H I) re- 
sembles the range of hydrogen column densities 
measured in QSO damped Lya systems (DLA 
hereinafter) - the systems with N(R I) > 2 x 10 20 
cm~ 2 , which are believed to originate in interven- 
ing galaxies or proto-galaxies located at cosmolog- 
ical distances (Wolfe et al. 1995). Thus we may 
expect to observe far-UV H 2 lines from electronic- 
vibrational-rotational bands shifted to the blue 
spectral range (A > 3000 A) at z > 2. The de- 
tection of the first H 2 -bearing cosmological cloud 
in the z arjs = 2.811 DLA system toward Q0528- 
250 was made by Levshakov & Varshalovich 
(1985) and later confirmed by Foltz, Chaffee, 
& Black (1988). Two additional H 2 systems 
were found at z a bs = 1-97 toward Q0013-004 
(Ge & Bechtold 1997) and at z abs = 2.34 toward 
Q1232+0815 (Ge & Bechtold 1999). The latter 
was recently confirmed by Petitjean, Srianand, 
& Ledoux (2000). The fourth identification of 
H 2 was made at z a b s = 3.39 toward Q0000-2620 
by Levshakov et al. (2000b, 2001) and a tentative 
fifth H 2 system was reported by Petitjean et al. 
(2000) at z abs = 2.374 toward Q0841+129. 

Molecular gas in the early Universe has also 
been detected through the spectral-line emission 
from the tracer carbon monoxide (CO) molecule 
in the ultraluminous galaxy IRAS 10215+4724 at 
z = 2.28 (Brown & Van den Bout 1992; Solomon, 
Downes, & Radfort 1992), in the clover- leaf quasar 
H1413+117 at z = 2.56 (Barvainis et al. 1997), in 
the quasar BR1202-0725 at z = 4.69 (Ohta et al. 
1996; Omont et al. 1996), in the quasar BRI1335- 
0418 at z = 4.41 (Guilloteau et al. 1997; Carilli, 
Menten, & Yun 1999), in the IR luminous galaxy 
J14011+0252 at z = 2.565 (Frayer et al. 1999), 
in the extremely red galaxy HR10 at z = 1.44 
(Andreani et al. 2000), in two powerful radio 
galaxies 4C60.07 and 6C1909+722 at z = 3.79 and 
3.53, respectively (Papadopoulos et al. 2000), in 
the lensed quasar APM08279+5255 at z = 3.911 
(Downes et al. 1999; Papadopoulos et al. 2001) 
and in the QSO SDSS 1044-0125 at z = 5.8 (Iwata 
et al. 2001). 

However, identification of H 2 lines in the 
damped Lya systems remains a rather rare 
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event. Most DLAs yielded only upper lim- 
its of /h 2 < 1(T 7 1Q- 6 . The lack of H 2 in 
these systems might be a consequence of the 
lower metal abundances, lower dust contents, 
higher temperature and ultraviolet fluxes when 
compared with the ISM in the Galaxy and in 
the Magellanic Clouds. Low molecular abun- 
dances found in high redshift DLAs resemble 
the properties of, at least, one local DLA sys- 
tem - the metal-deficient (Z/Zq ~ 1/50) star- 
burst blue compact galaxy I Zw 18 with N(R I) 
~ 3 x 10 21 cm~ 2 (Kunth et al. 1994) where re- 
cent FUSE observations set a limit of /h 2 -C 10 -6 
(Vidal-Madjar et al. 2000). 

The measured upper limits of /h 2 in high red- 
shift metal-poor DLAs are slightly lower than the 
calculated fraction of molecules in the pristine gas 
- /h 2 — 10 -6 -j- 10~ 5 depending on cosmological 
models (e.g. Galli & Palla 1998). Although UV 
photons produced by the first objects can quickly 
destroy all H 2 in the pristine gas and thus de- 
lay further structure formation (Haiman, Rees, & 
Loeb 1997), we may expect that the electrons in 
the photoionized regions around first stars will in 
turn act as a catalyst in the H~ channel to pro- 
duce H 2 at redshifts z < 50 where the H~ chan- 
nel proceeds much more quickly than the H + one 
(Fuller & Couchman 2000). Therefore direct mea- 
surements of gas density, UV radiation field, and 
H 2 formation and destruction rates in DLA sys- 
tems may shed light on star formation processes 
at very high redshift. 

In this paper we present results of the study 
of H 2 lines identified in the sixth DLA system at 
z abs = 3.025 toward Q0347-3819. We also use the 
information from the molecular lines together with 
H and metal lines to refine the deuterium abun- 
dance determination by D'Odorico, Dcssaugcs- 
Zavadsky, & Molaro (2001, hereinafter DDM) and 
the values of metal abundances (DDM; Prochaska 
& Wolfe 1999). 

The outline of the paper is as follows. We dis- 
cuss observations in Section 2. The methods of 
data analysis are presented in Section 3. The iden- 
tification of H 2 lines and the calculation of their 
proper redshift is given in §3.1. The broadening 
^-parameter is evaluated in § 3.2. The column den- 
sities N(3) for the rotational levels from J — to 5 
of the lowest vibrational level v = in the ground 
electronic state (for the H 2 notation, see 



Field et al. 1966), the fractional molecular abun- 
dance, and the ortho:para-H 2 ratios are calculated 
in § 3.3. In § 3.4 we compute the excitation tem- 
perature T cx . The results of a stacking technique 
applied to the measurements of the abundances 
of Zn II and Cr II ions and the dust-to-gas ratio 
in the whole z^bs = 3.025 system are estimated in 
§3.5. In § 3.6 we analyze within one self-consistent 
model all available profiles of H, D, neutral atoms 
and low ions to infer column densities and relative 
abundances of D/H and metals. In Section 4 we 
discuss the results and present their further anal- 
ysis. In §4.1 we estimate the kinetic temperature, 
Tkin, and the gas volumetric density, nn, using 
the C ir/C II and C I/C II ratios. The distri- 
bution of H 2 over the rotational levels, the photo- 
dissociation rate of H 2 and the rate of the molec- 
ular hydrogen formation upon grain surfaces are 
considered in §4.2. The revealed [a-element /iron- 
peak] enhancement is discussed in § 4.3. In Section 

5 we summarize the results. 

2. Observations 

Q0347-3819 is a bright quasar (V = 17.3) with 
z cm = 3.23 discovered by Osmer & Smith (1980), 
which shows a damped Lya system at z a b s = 
3.0245. A column density of N(H I) = (6.3±1;|) x 
10 20 cm -2 was first measured by Williger et al. 
(1989) and later confirmed by Pettini et al. (1994) 
as N(H I) = (5 ± 1) x 10 20 cm" 2 . Q0347-3919 
has a rich absorption-line spectrum studied by 
Pierre, Shaver, & Robertson (1990), Centurion et 
al. (1998), Ledouxet al. (1998), and by Prochaska 

6 Wolfe (1999). 

High-quality data of this QSO in the UV (3650 
A< A < 4900 A) and in the near-IR (6700 
A < A < 10000 A) ranges were obtained during 
the commissioning of the UVES on the VLT 8.2m 
Kueyen telescope at Paranal, Chile, in December 
1999, and have been released for public use. The 
UVES spectrograph is described by D'Odorico et 
al. (2000) and details of the data reduction arc 
given in Ballcster et al. (2000). Two exposures 
of 4500 s and 5000 s, covering the UV and the 
near-IR ranges, were obtained with a resolution 
FWHM ~ 7.0 km s _1 and ~ 5.7 km s _1 , respec- 
tively. An average signal-to-noise ratio per pixel 
of ~ 20, 25, 40 and 15 was achieved in the final 
spectrum at A ~ 3650 A , 4900 A , 6700 A , and 
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10000 A , respectively. 
The 

^abs — 3.025 DLA exhibits a multicompo- 
nent velocity structure with the dominating com- 
ponent at z a b s = 3.024856 as found by DDM who 
discussed the first detection of the deuterium lines 
in DLAs. The total H I column density de- 
rived in DDM from the best fit of Ly/3 equals 
(4.3 ± 1.0) x 10 20 cm~ 2 , which agrees well with the 
measurement of Pettini et al. (1994) from Lya. 
The ionization structure of this DLA is also com- 
plex. There are at least two gas components: a 
warm gas seen in lines of neutral atoms, H 2 , and 
low ions, and a hot gas where the resonance dou- 
blets of C IV and Si IV are formed. These high- 
ion transitions, observed by Prochaska & Wolfe 
(1999) with the Keck HIRES spectrograph show, 
however, different profile shapes as compared with 
the low ion-transitions. The physical properties of 
the hot gas component will not be studied in the 
present paper. 

3. Measurements 

All absorption lines studied in the present pa- 
per lie in the Lya forest where their analysis is 
complicated by possible H I contaminations and 
uncertainties in a local continuum, drawn by us- 
ing 'continuum windows' in the Lya forest which 
could be located quite far from the feature under 
analysis. This kind of uncertainties is illustrated in 
Figure 1 where a portion of the Q0347-3819 spec- 
trum with a few identified H2 lines is shown. To 
achieve the highest possible accuracy under these 
circumstances, we have always used in parallel sev- 
eral different lines to derive the physical parame- 
ters of the DLA system. 

3.1. Redshift of the H2 system 

Wavelengths, oscillator strengths, and radia- 
tion damping constants for the Lyman and Werner 
bands of the H2 molecule were taken from the the- 
oretical calculations of Abgrall & Roueff (1989). 
These data seem to be reasonably accurate as 
shown by high spectral resolution IMAPS mea- 
surements of molecular hydrogen toward the Orion 
belt stars (Jenkins et al. 2000). 

The identification of the H2 lines was started 
from the subsystem which has the largest H I col- 
umn density (see § 3.6). Preliminary inspection of 
the Q0347-3819 spectrum revealed 88 absorption 



features at the expected positions of the H2 lines. 
However, not all of them were suitable for further 
analysis due to H I contaminations. For this rea- 
son, to estimate z-r 2 we used only those H2 lines 
which show symmetric profiles with their blue and 
red sides reaching the continuum. These lines are 
listed in Table 1 (see also Figures 5 and 6). 

Given a fit to the continuum and the line mea- 
surement interval, the equivalent width W\, the 
line center (A) and their errors <Jw x arL d cr (A) can 
be computed from the following equations (e.g. 
Levshakov et al. 1992): 
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where A A is the wavelength interval between pix- 
els, M is the number of pixels involved in the 
sum, Fi is the normalized intensity in the ith pixel 
whose wavelength is Xi, (Ji is the noise level in 
the ith pixel, W = Wa/AA is the dimensionless 
equivalent width in pixels, and (cr c /C) is the mean 
accuracy of the continuum fit over the line mea- 
surement interval. 

The errors in W\ and (A) are caused mainly by 
the finite number of photons counted. In addition, 
for weak absorption features (W <C 1) the second 
term in equation (3) is of about the same order 
of magnitude as the first one, meaning that the 
identification of weak H 2 lines is very sensitive to 
the accuracy of the local continuum fitting. 

The analyzed 35 transitions of the Lyman and 
Werner bands arc listed in column (2) of Ta- 
ble 1, whereas their J-values are given in col- 
umn (1). The measured heliocentric position of 
the center of each spectral feature and its stan- 
dard deviation were calculated according to equa- 
tions (2) and (4) and are presented in columns 
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(3) and (4), respectively. The observed equiva- 
lent widths (column 5) and their errors (column 
6) were calculated according to equations (1) and 
(3). Note that the rest frame equivalent width 
equals WA,rost = W\/{\ + Zn 2 )- In these calcula- 
tions we set for all lines (a c /C) = 1% to account 
for the uncertainty in the local continuum level. 
A formal mean value of (a c /C) estimated from 
the continuum windows is less than or about 1%. 
In particular, (a c /C) = 0.70% for L3-0R(0) and 
1.03% for L7-0P(5) which are shown in Figures 5 
and 6, respectively. Column (7) lists the measured 
redshifts for each H2 transition, their errors | a z \ 
= I <T(x) I /Ao (where Ao is the rest frame wave- 
length from Abgrall & Roueff 1989) are shown in 
column (8). The weighted mean rcdshift yields 
z H2 = 3.024855 ± 0.000005, which coincides ex- 
actly with the redshift of the strongest main com- 
ponent in metal line profiles given by DDM. Thus, 
at a given S /N ratio and spectral resolution we do 
not find any systematic shifts larger than a few 
km s _1 in the radial velocities between neutral 
atoms and H2 molecules in the molecular cloud 
at z a b s = 3.025. 

Previous Copernicus investigations of diffuse 
gas in the Galaxy showed that the mean radial 
velocity of H 2 lines differs from the correspond- 
ing mean of the atomic lines. For instance, the 
differences of Av ~ 5 km s _1 toward A Sco 
(Spitzer & Jenkins 1975), Aw ~ 11 km s _1 and 
~ 6 km s _1 toward, respectively, r Sco and ( Oph 
(Spitzer, Cochran, & Hirshfcld 1974) were de- 
tected. At the UVES spectral resolution such dif- 
ferences - in case they were present - could be 
easily distinguished. 

3.2. Line broadening 

The measurement of the broadening b 7 param- 
eter and the calculation of the confidence ranges 
for its mean value were carried out in two steps. 

First, a least-squares analysis was applied to 
the line fitting. We combined 17 H2 lines (L3-0R0, 
L7-0R0, W3-0Q1, L3-0P1, W0-0R2, W0-0Q2, L3- 
0R2, L4-0R2, L2-0P3, L2-0R3, L4-0R3, L4-0P3, 
L8-0P3, L12-0R3, L4-0P4, L8-0P4, W1-0Q5) from 
Table 1 and fitted them simultaneously with the 



By definition b = y2a v , where u v is the one-dimensional 
Gaussian velocity dispersion of molecules along the line of 
sight. 



individual Gaussian profiles. These lines were cho- 
sen because they exhibit 'clean', uncontaminated 
profiles. The H 2 line centers were fixed at the 
mean zh 2 and the theoretical profiles were con- 
volved with the instrumental point-spread func- 
tion (PSF) which was approximated by a Gaus- 
sian with FWHM = 7.0 km s _1 , i.e. 6 PSF = 4.2 
km s _1 . The Ax 2 method (Press et al. 1992) was 
used to estimate the most probable b value and its 
error 8 . For a given current value of b taken from 
the range [& mm , &max], we minimized \ 2 (normal- 
ized per degree of freedom) by varying the line in- 
tensities independently (we consider Xmin as a X 2 
distribution with v = M — M degrees of freedom, 
where the number of data points N = 192 and the 
number of fitted parameters M. — 6 for this par- 
ticular case). With this we tried to check possible 
H I contaminations : the optimization procedure 
would fail to reach Xmin — 1 if the chosen H2 lines 
were essentially contaminated by the Lya forest 
absorption. The calculated Xmin values in depen- 
dence of b and the la and 2a confidence levels are 
shown in Figure 2. With xLn = °- 925 the most 
likely value for 6h 2 was found to be 2.80 ± 0.45 
km s _1 (s.d.). 

In the second step, to check the obtained value 
and to infer the Doppler parameter independently, 
we calculated an optimal curve of growth for the 
H2 lines from Table 1. The curve is shown by 
the solid line in Figure 3. In principle, it can- 
not be excluded that some unknown absorption 
lines blend with H 2 features, despite their sym- 
metric and undisturbed profiles. Thus, the scatter 
of points in Figure 3 can be explained by this ef- 
fect as well as by errors remaining after the contin- 
uum fitting. It is also seen from Figure 3 that the 
strongest absorption lines from J = 1 and 3 levels 
are slightly saturated and therefore their equiv- 
alent widths are sensitive to the assumed broad- 
ening model. The data for these levels give the 
effective Doppler width b ~ 3 km s _1 . We also de- 
picted two curves for b = 2 km s _1 and 4 km s _1 to 
demonstrate the complete consistency of 6-values 
estimated in both approaches. 



In the present paper, the estimated standard errors are 
slightly larger than those calculated by the Monte Carlo 
method because the pipeline resampling to constant wave- 
length bins introduces a correlation between data point val- 
ues, whereas the estimated mean values of the physical pa- 
rameters remain unbiased (see Appendix A). 
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We note that Jenkins & Pcimbcrt (1997) have 
measured an interesting phenomenon of increas- 
ing 6-values for H2 lines arising from progres- 
sively higher rotational levels in a spectrum of 
( Ori A. The most probable Doppler parameters 
were found to be : 6(0) = 0.8 km s" 1 , 6(1) = 2.7 
km s" 1 , 6(2) = 5.3 km s" 1 , 6(3) = 5.9 km s"\ 
and 6(5) = 9.0 km s _1 . As they noted, such profile 
shapes are inconsistent with a simple model of UV 
optical pumping which governs the H2 rotational 
populations in an optically thin, homogeneous 
medium. They suggested that these H 2 profiles 
with changing shapes 'may represent molecules 
forming in the postshock gas that is undergoing 
further compression as it recombines and cools'. A 
similar phenomenon may also be observed in high 
rcdshift DLAs which can be associated to star 
forming regions where shocks are very probable. 

The z a b s = 3.025 system presents the weakest 
members of the Lyman and Werner series with 
J = 4 and 5 but, unfortunately, we arc unable 
to compare their 6-values with those from lower 
rotational transitions because our signal-to-noise 
ratio is not sufficiently high. Thus, in the follow- 
ing we assume a simple one-component model for 
the molecular cloud with a common broadening 
parameter 6h 2 = 2.80 km s _1 for all molecular 
hydrogen lines. 

3.3. H 2 abundances 

In the final fitting of each line profile we com- 
bined together the H 2 lines with the same J to 
yield the corresponding N(J) column densities and 
their confidence regions. A x 2 minimization pro- 
cedure was applied simultaneously to all lines from 
a given rotational level J. To control solutions at 
this stage, we added to the lines listed in Table 1 
a few weaker transitions where we observe con- 
tinuum without any pronounced absorption fea- 
tures at their expected positions. The calculated 
Xmin values and the la and 2a confidence levels 
for each individual rotational level J are shown 
in Figure 4. The curves x 2 [iV(J)] were computed 
with the fixed mean 6h 2 = 2.80 km s _1 and the 
mean zh 2 . The H 2 lines involved in the analysis 
are as follows: L3-0R0, and L7-0R0 (TV = 22); 
L3-0P1, and L14-0P1 (Af = 20); W0-0R2, W0- 
0Q2, L3-0R2, L4-0R2, and L5-0R2 (Af = 51); W2- 
0R3, L2-0R3, L2-0P3, L3-0R3, L3-0P3, L4-0R3, 
L6-0R3, L7-0P3, L8-0P3, L12-0R3, and L15-0R3 



(Af = 83); L5-0R4, L6-0R4, and L8-0P4 (Af = 34); 
L7-0P5, and L12-0R5 (Af = 21). The most prob- 
able column densities and their errors at the la- 
confidence level are: N(0) = (1.90 ± 0.35) x 10 13 
cm" 2 , N(l) = (1.60 ±0.18) x 10 14 cm" 2 , N(2) = 
(5.40 ±0.35) x 10 13 cm" 2 , N(3) = (1.22 ±0.05) x 
10 14 cm" 2 , A(4) = (2.60 ± 0.43) x 10 13 cm" 2 , 
and N(5) = (2.90 ± 0.70) x 10 13 cm" 2 . From 
that we calculate the total column density in H 2 , 
JV(H 2 ) = £j N(J) = (4.10 ±0.21) x 10 14 cm" 2 , 
and the ratio of the total column densities in the 
ortho- and para-H 2 states N OIt ho/N paia — 3.1±0.3 
(s.d.). This result can be compared with two other 
molecular clouds at z abs = 2.8112 (Q0528-250, 
Srianand & Petitjean 1998) and z abs = 2.33771 
(Q1232+0815, Srianand et al. 2001) where N(3) 
values were also measured for different J. Using the 
published A^(J) values and their errors, we found 
AUho/^Vpara = 3.1 ±0.8 (s.d.) at z abs - 2.8112 
and JVortho/iVpara = 2.4 ± 0.9 (s.d.) at z abs = 
2.33771. 

According to our calculations presented in Sec- 
tion 3.6, the total H I column density in the z abs — 
3.025 system is (4.225 ± 0.045) x 10 20 cm- 2 . With 
iV(H 2 ) = (4.10 ± 0.21) x 10 14 cm" 2 , the ratio of H 
nuclei in molecules to the total H nuclei is 

/h2 = ^^hT = (l94±0 - 10)x10 "- 

To test the consistency of the estimated physi- 
cal parameters with the observed absorption fea- 
tures, we calculated theoretical profiles for all 88 
H 2 lines using the derived mean zh 2 , 6h 2 , and 
N(J) values. Besides we calculated 8 profiles for 
H 2 transitions where no absorption features have 
been observed at the continuum level. The latter 
are L3-0R4, L9-0R4, L4-0R5, L9-0R5, L11-0P5, 
W0-0Q5, L6-0P6, and L7-0R7. The last two of 
them give upper limits to A^(6) < 2.0 x 10 13 cm -2 
(la c.L), and to N(7) < 8.0 x 10 12 cm" 2 (la c.L). 

The observational and over-plotted theoretical 
H 2 profiles are shown in Figures 5 and 6. In these 
figures possible blends with other lines from the 
^abs = 3.025 DLA and from the metal systems 
listed in Appendix B are also indicated. We found 
only one large inconsistency between observable 
and theoretical profiles, namely for the W2-0Q(3) 
line (Figure 6). But in this case a very wide 
absorption trough between 3901 A and 3921 A 
redward the redshiftcd position of this molecular 
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line (A bs = 3900.3 A) prevents the accurate es- 
timation of the local continuum level. Uncertain- 
ties in the local continuum level calculations may 
also cause over-absorption in the synthetic profiles 
as compared with the observable profiles seen in 
Wf-ORO+Rf, W2-0Q1 (Figure 5), W0-0Q3, L3- 
0P3, L4-0P3, L10-0R3, L12-0R3, W0-0Q5, and 
L4-0R5 (Figure 6). In all other cases the synthetic 
H 2 profiles are in agreement with the absorption- 
line features identified with the H2 lines in the UV 
spectrum of Q0347-3819. 

We have also tried to detect the H 2 molecules 
in the blue component at Av ~ — 17 km s _1 (see 
§ 3.6) with respect to the main H 2 -bearing com- 
ponent but we found only an upper limit. The ap- 
plied method is similar to that employed by Lev- 
shakov et al. (1989) to improve the detectability of 
weak carbon monoxide (CO) UV lines in damped 
Lya systems. 

The procedure consists of shifting the observed 
spectral regions containing expected positions of 
the H 2 transitions, which are not contaminated by 
the apparent Lya forest absorption, into the rela- 
tive velocity scale centered at z — 3.024637, - the 
rcdshift of the blue subcomponent of the dominant 
metal absorption in the z a b s = 3.025 DLA, and of 
the consequent stacking of these regions. Before 
co-adding each section, the data are rebinned with 
a bin size of 1.45 km s _1 . The spectra were then 
averaged with weights in accord with the squares 
of the expected relative line strengths divided by 
the respective signal-to-noise ratios. 

The result obtained from this stacking proce- 
dure is shown in Figure 7. The composite spec- 
trum formed by the addition of 18 H 2 bands (with 
no pronounced absorption features in the range 
Av = ±6 km s _1 ) which are labeled in the corre- 
sponding panels is depicted in the bottom panel. 
The resulting co-added data have an average S/N 
~ 100 and show no absorption at v = km s _1 , 
whereas the blue wing of the composed molecular 
line from the z = 3.024855 system is clearly seen 
at v ~ +16 km s _1 . 

The composite spectrum gives Merest < 0.6 
mA (3crc.l.) for the equivalent width of the sum 
of the H 2 bands in the rest frame of the absorber. 
This limit is sensitive, in general, to errors in con- 
tinuum fitting. In our case the resulting profile be- 
tween —6 km s _1 < v < 6 km s _1 has no apparent 
deviations from unity, the estimated mean value 



is (C) = 1.003 ± 0.003 and thus (a c /C) = 0.003. 
Substituting these values in equation (3) we obtain 
the limiting equivalent width (W <C 1) <y\im = 0.2 
mA and TUA, rcs t = 0.6 mA (3(7i im ). 

For a weak absorption line with the central op- 
tical depth r <C 1 , the column density is indepen- 
dent of the broadening parameter b and is linearly 
proportional to the equivalent width W\: 

7V = 1.13xl0 20 -^,cm- 2 (5) 

where / is the oscillator strength for the absorp- 
tion transition, Ao is the line center rest wave- 
length in A, and W\ is the rest equivalent width 
in A. 

Using this linear approximation to the curve 
of growth, we calculate the composite equivalent 
width Wj for the H 2 bands from the same J level 

Wj = aN 3 J2^,ihi, (6) 

i 

where a = 8.85 x 10 -21 . 

Then the total equivalent width for all H 2 bands 
presented in Figure 7 is 

We = W * = a Yl ej Nj = a ' ( ? ) 
j j 

where ej = J2i ^ i s an effective value (in 

our estimation e = maxej), and JVh 2 = J2j ^J- 
The coefficients ej do not vary significantly and 
for the H 2 bands from Figure 7 e = 6.1 x 10 4 , 
ei = 5.5 x 10 4 , e 2 = 5.3 x 10 4 , and e 3 = 9.8 x 
10 4 . Thus, the aforementioned upper limit W\ < 
0.6 mA and equation (7) yield N(R 2 ) < 7.0 x 
10 11 cm" 2 (3cr c.l.) and / Ha < 3.3 x 10~ 9 (3cr 
c.L), which is well below that found in Galactic 
H I clouds. It means that molecular hydrogen is 
not homogeneously distributed in the z a b s = 3.025 
absorbing region. 

The measured iV or tho/^para ratios in the DLAs 
at z abs = 3.025 (Q0347-3819), 2.8112 (Q0528- 
250), and 2.33771 (Q1232+0815) are equal to 3 
(within the uncertainty interval), which is just 
the equilibrium value for the ortho:para ratio at 
the hot formation of H 2 on grain surfaces in rel- 
atively transparent clouds (e.g., Spitzer & Mor- 
ton 1976). In the pristine gas, the ortho:para 
H 2 ratio decreases from the equilibrium value of 
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3 at z ~ 100 to its 'freeze-out' quantity of 0.25 at 
z < 20 (Flower & Pineau des Forets 2000). Thus, 
we can conclude that H2 observed in these DLAs 
has undergone a re-forming phase following photo- 
dissociation of the pristine H 2 . 

3.4. Excitation temperature 

In the ground clectronic-vibrational state, all 
the rotational levels with odd J are nuclear triplet 
states ( ortho-H 2 ), and all those with the even 
J are singlet states ( para-H2 ). Thus, for the 
ortho-H 2 the statistical weight, g(3), of a level J 
equals 3(2J+1), whereas for the para-molecules 
<?(J) = (2J + 1). In thermal equilibrium, the colli- 
sionally dominated distribution of molecules over 
different J levels is given by a Boltzmann law with 
an effective excitation temperature T cx 



N(3) _ g(J) 
N(0) g(0) 



exp 



B v 3(3 + 1) 



(8) 



where the constant B v equals 85.36 K for the vi- 
brational ground state. 

This equation shows that T ex is proportional 
to the negative inverse of the slope of a line 
(called 'the excitation diagram') drawn through 
the points of the respective J levels in a plot 
log [N(3)/g(3)} versus E(3), where iV(J) is the 
column density measured in the rotational level 
J, E(3) is the excitation energy of this level rel- 
ative to J = 0, and g(3) is its statistical weight. 
Figure 8 shows this diagram for different J lev- 
els. A least-squares fit, weighted by 1/er 2 for each 
log [N(3)/g(3)\, yields T cx = 825 ± 110 K (s.d.) 
and log N(0) = 13.31 ± 0.03 in cm" 2 (s.d.). This 
value of N(0) agrees well with what has been mea- 
sured from profile fitting in § 3.3. Thus, we find 
that the population of the six lowest J levels can 
be fitted with a single Boltzmann distribution. 

The unique T cx for all levels from J = to 5 may 
imply that the H2 rotational population is domi- 
nated by collisions in a gas whose kinetic temper- 
ature is about 800 K. The kinetic temperature in 
the H2 cloud can be estimated from the widths 
of the H2 and metal lines. In § 3.6, we show that 
profiles of all metals from C II to Zn II can be 
described with the same broadening parameter of 
b Mc = 2.71 ±0.05 km s" 1 . With 6 H2 = 2.80 ±0.45 
km s _1 , it gives an upper limit of Tkin < 430 K 
in the main H2-bearing component. The fact that 



Ikin is less than T cx indicates that the population 
of the low J levels is caused not only by collisions 
but by UV-pumping as well. The J = and 1 
rotational levels, being usually populated in dense 
molecular clouds by thermal proton collisions, may 
not be in equilibrium with the kinetic tempera- 
ture in lower density clouds which are not optically 
thick in the L0-0R(0), R(l), P(l) and W0-0R(0), 
R(l), Q(l) lines. In low density clouds the self- 
shielding effects are negligible and do not reduce 
significantly UV-pumping in these lines (e.g., Shull 
& Beckwith 1982). For the further analysis we set 
Tidn = 400 K to be consistent with the measured 
b- values. 

3.5. Cr/Zn and dust content 

Copernicus H2 survey has shown that in our 
Galaxy the molecular fraction /h 2 is correlated 
with E(B — V) and with the total neutral hydro- 
gen column density (e.g. Spitzer & Jenkins 1975). 
The presence of dust grains in damped Lya sys- 
tems is usually estimated from the abundance ra- 
tio [Cr/Zn] 9 assuming that Zn is undepleted (Pet- 
tini ct al. 1994). 

Preliminary visual inspection of the red part 
of the Q0347-3819 spectrum has shown weak ab- 
sorption features at the expected positions of Zn II 
and Cr II lines. In order to increase sensitivity 
and to measure the abundance ratio [Cr/Zn] we 
applied a stacking procedure similar to the one 
described in §3.3 to the Zn II A2026.136 A and 
A2062.664 A lines, and to the Cr II A2056.254, A 
A2062.234, A and A2066.161 A lines. The wave- 
lengths and the oscillator strengths for these lines 
are taken from Morton (1991), and from Berge- 
son & Lawler (1993), respectively. The resulting 
composite spectra of Zn II and Cr II were in- 
cluded in the metal analysis in the next § 3.6 (see 
Figure 10 and Tables 3 and 4). In this Figure, 
two pronounced absorption features at Av ~ +6 
km s- 1 (Cr II) and at Av ~ +13 km s" 1 (Zn II) 
are telluric lines. The former partly blends the 
red wing of the Cr II line, therefore we do not use 
this portion of the Cr II profile in the fitting pro- 
cedure. The estimated total column densities are 



9 [X/Y] = (X/Y) - (X/Y) , where (X/Y) is the logarithmic 
value of the element ratio by number without reference to 
the solar value. Solar abundances (X/Y)q are taken from 
Grevessc & Sauval (1998). 
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7V(Zn II) = (6.19±0.51) x 10 11 cm- 2 , and iV(Cr II) 
= (2.76 ±0.22) x 10 12 cm" 2 . With these numbers 
one finds [Cr/Zn] = -0.37 ± 0.06. 

In Table 2 we summarize the results of the 
abundance measurements of molecular hydrogen, 
[Cr/Zn] (and/or [Fc/Zn]), and neutral hydrogen 
in the DLAs available up to now. The pres- 
ence of dust in these DLAs and the relation be- 
tween the fractional molecular abundances and 
relative heavy clement depletion is illustrated in 
Figure 9. The measured data for the five ex- 
tragalactic molecular clouds (marked by numbers 
in Figure 9) can be fitted with a linear law : 
log/ H2 = A [Cr/Zn] + 5, with A = -6.053 ±0.116 
and B = —7.888 ± 0.051. The linear regression is 
shown in Figure 9 by the solid line 10 . Again, simi- 
larly to our previously obtained results (Levshakov 
et al. 2000b), we find a strong correlation be- 
tween log /h 2 in the DLA systems with detectable 
H 2 absorption and the heavy element depletion. 
The presence of upper limits below the linear re- 
gression line may be explained by the suggestion 
that in these systems a relatively high temperature 
of dust grains prevents molecular formation. 

The degree of the heavy element depletion may 
be used to estimate the dust-to-gas ratio, k, if two 
measured elements have the same nucleosynthesis 
history (cf. Vladilo 1998) : 



k = 



10[ Zn / H lob= 



/zn.ISM — /Cr,ISM 



(l [Cr/Zn] obs _ ^ ) (9) 



where the fraction in dust /z n ,iSM = 0.35 and 
/cr.iSM = 0.92 refers to the Galactic interstellar 
medium. 

Error propagation method applied to the 
column density measurements yields [Zn/H] = 
-1.503 ± 0.054, and then Eq. (9) results in 
the average value for the whole DLA /cdla = 
0.032 ±0.005 (la c.l.) showing that the dust con- 
tent in the z a b s = 3.025 absorber is approximately 



10 For the published 27 DLAs where all three elements Zn II, 
Cr II, and Fe II were measured, [Cr/Fc] = 0.12 ± 0.13 
(s.d.). This implies that Fe is slightly more depleted than 
Cr assuming both elements track each other in the abun- 
dances and, thus, the slope in our regression analysis may 
be slightly steeper. However, the accuracy of the [Cr/Zn] 
and [Fc/Zn] measurements is not sufficiently high to pro- 
vide an unambiguous rcnormalization of one element to the 
other. Therefore [Fc/Zn] and [Cr/Zn] are used interchange- 
ably in the present calculations. 



30 times lower as compared with the Milky Way. 
If, however, we consider the same quantity for the 
H 2 -bearing cloud using the data from Table 3, 
then k ~ 0.047, i.e. the dust content in the molec- 
ular cloud is about 1.5 times higher as compared 
with the average value. 

3.6. D/H and metal abundances 

In this section we study metal profiles and de- 
rive a new D/H ratio by fitting simultaneously 
the available Lyman series lines and metals. The 
knowledge of the proper distribution of the veloc- 
ity components along the line of sight is especially 
crucial for the D/H measurements, since D I lines 
are always blended with the blue wings of the cor- 
responding H I lines (for numerous examples, see 
Levshakov, Kegel, & Takahara 1999 and references 
cited therein). 

In the UVES spectrum of Q 0347-3819, at 
Zubs = 3.025 we have identified absorption lines 
of C II, N I, O I, Si II, P II, Ar I, Cr II, Fe II, 
Zn II, and C II*. The identified metal absorption 
lines are shown in Figure 10, and in Figure li- 
the atomic hydrogen lines. The atomic data were 
taken mainly from Morton (1991) with a few addi- 
tions of new results for Si II from Spitzer & Fitz- 
patrick (1993) and Charro & Martin (2000), for 
Ar I from Federman et al. (1992), for Cr II and 
Zn II from Bergeson & Lawler (1993), for Fe II 
from Cardelli & Savage (1995), Raassen & Uylings 
(1998), and Howk et al. (2000). 

To construct a model for the radial velocity dis- 
tribution, we firstly analyzed the profiles of each 
individual atom/ion to select lines (or parts of 
them) which are less contaminated and disturbed 
by the Lya forest absorption. The chosen profiles 
and/or their components are marked in Figures 10 
and 11 by horizontal lines at the bottoms of the 
panels (these horizontal lines also show the pix- 
els involved in the final optimization procedure). 
Other spectra in Figures 10 and 11, not marked 
by the horizontal lines, were used to control con- 
sistency of the obtained results. 

Our model is based on the following assump- 
tions: (i) the whole DLA is considered as a contin- 
uous absorbing region with fluctuating gas density 
and velocity fields, (ii) The metallicity for each el- 
ement is set constant through the DLA. (Hi) Some 
of the observed species may have similar profiles, 
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i.e. their fractional ionization ratios remain con- 
stant along the sightlinc. In this case complex 
profiles can be modeled using the sum of the same 
number of components with identical centers and 
equal relative intensities for any pair of lines (Lev- 
shakov et al. 1999b). 

The line profiles are described by the sum of m 
Voigt functions : 

m 

ta = E N i V K A - A o - Ao-)/AAi] , (10) 
»=i c 

where t\ is the optical depth at wavelength A 
within the line profile, Ni is the column density 
in the ith component, Vi is its Dopplcr shift, and 
AAi = Ao&i/c is its Dopplcr width with c being 
the speed of light and Ao the line center. 

Equation (10) requires additional suggestion 
about the broadening mechanism. We assume 
that the main broadening is caused by bulk mo- 
tion. For the component seen in the H2 lines this 
suggestion is supported by the measured equiva- 
lence of 6jj 2 and &Me> whereas for other compo- 
nents it is a model assumption which has to be 
checked during the optimization procedure. 

Now we can describe our model in detail. It is 
fully defined by specifying {Ni}f =1 - the column 
density of the 1st components for each ion, with k 
being the total number of ions ; {Tlu}™ =2 ~ the rel- 
ative intensities of the subcomponents, where for 
a given ion TZu = Ni/N\ ; {wi}™! - the Dopplcr 
shifts ; and 1 - the Dopplcr widths. All 

these parameters are components of the parame- 
ter vector 9 = {61,62, . . . ,6 q }. In addition to 6, 
the D/H ratio is to be set. To estimate 9 from the 
absorption line profiles, we minimize the objective 
function 

L M e 

x 2 (d/h) = - E E > 

1=1 1=1 

(ii) 

where D/H is the total hydrogen isotopic ratio for 
the whole system. In (11), Tff. is the observed 
normalized intensity of the spectral line I, and cr^.i 
is the experimental error within the zth pixel of the 
line profile. !FjF\. (6) is the calculated intensity of 
line I at the same ith pixel having wavelength A, 

^,1 (0) = J e~ T ^ 0p SF (A, - A) dX , (12) 



where t^a^) is defined through Eq.(10). The in- 
strumental point-spread function </>psf(A) is ap- 
proximated by a Gaussian with FWHM = 7.0 
km s- 1 (see § 3.2). 

The total number of spectral lines involved in the 
optimization procedure is labeled by L, and the 
total number of data points M — J2e Me, where 
Mi is the number of data points for the £th line. 
The number of degrees of freedom is labeled by v. 

After preliminary trial computations we were 
able to separate all available lines into two groups 
consisting of similar profiles (the fitting of each 
group of lines yielded acceptable x 2 values) . Sim- 
ilarity in this context means that the density 
weighted radial velocity distribution functions 
for elements combined in a group are equal [see 
eq.(15) and explanations below]. First group in- 
cludes k = 9 different ions (H I, C II, C II*, N I, 
O I, Si II, P II, Ar I, and Zn II) with the total 
number of data points M. = 566 and v = 544 (we 
consider C II* as a separate element since it has 
a column density different from C II). The second 
one consists of k = 2 ions (Fe II and Cr II) with 
M = 110, v = 99. These two groups are essen- 
tially the refractory and non- (or little) refractory 
elements. 

We do not consider the absorption features seen 
at Av ~ 40 and 60 km s" 1 in the C II A1036.3 A , 
N I A1134.2 A, O I A1039.2 A, and Si II A1193.3 
A profiles (Figure 10) since we are mostly inter- 
ested in the main H 2 -bearing component (Av = 
km s _1 ) and in the blue side parts of the pro- 
files which restrict the configuration of the veloc- 
ity field important for the analysis of the observed 
D I lines. For this reason only blue parts of the 
hydrogen H I Lya, Ly-5, Ly-8, Ly-10, and Ly-12 
lines were included in the analysis (see Figures 11 
and 12). The components at Aw ~ 40 and 60 
km s _1 increase the total column densities mea- 
sured in the range —50 km s _1 < v < 80 km s _1 
at about 8% only. 

The estimation of the parameter vector 9 by 
means of the x 2 -minimization is a typical ill-posed 
problem (Tikhonov & Arsenin 1977) because of 
the presence of two smoothing operators in (10) 
and (12) which prevent from deducing a unique 
value of 9. To stabilize the solution, we augment 
X 2 by a regularization term in a way similar to that 
discussed in Levshakov ct al. (1999b) - instead of 
minimizing % 2 , defined in (11), we now seek the 
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minimum of another objective function, 

£ a (6)= X 2 (0)+am, (13) 

where ip is a penalty function and a the regular- 
ization parameter. 

The choice of t/j for nonlinear problems oc- 
curs rather heuristically. For our particular DLA 
we used additional information obtained from the 
analysis of H2 . Namely, we searched for the min- 
imum of C a (6) in the vicinity of the main compo- 
nent (i.e. Ai>h 2 — km s _1 ) which is supposed to 
have b ~ &h 2 - Then the easiest choice for ip will 
be 

= (— —) +(Av-Av H2 ) 2 , (14) 

and we set the regularization parameter a = 1. 

The importance of the regularization is illus- 
trated in Figure 13, where different velocity distri- 
bution functions p(v) are shown. The differences 
in shapes of p(v) may cause different D/H val- 
ues. For instance, being analyzed separately, the 
H+D profiles yield a shallow and smoother p(v) 
(the dashed curve in Figure 13), whereas a simul- 
taneous fitting with metals reveal a sharp spike at 
Av ~ km s _1 (the solid curve) and a slightly dif- 
ferent structure of the blue wing of p(v) at Av < 
—5 km s _1 . The difference in D/H is, indeed, es- 
sential: D/H = (2.24±0.67) x 1(T 5 in the former 
case (see DDM), and D/H = (3.75 ±0.25) x 1CT 5 
in the latter one (see Figure 14). 

The estimated column densities together with 
their errors and relative abundances in the H 2 - 
bearing cloud are listed in Table 3. The errors 
were calculated by using the A% 2 method. They 
reflect internal uncertainties of the physical pa- 
rameters within the adopted model. Total col- 
umn densities (calculated in the range —50 km s _1 

< v < 20 km s -1 ) for all elements, their relative 
abundances and corresponding errors are listed in 
Table 4. Besides, we also estimated in the main 
component la upper limits to N(C I) < 4.2 x 10 11 
cm" 2 from the C IA1277.2 A line, to N(0 I*) 

< 8.2 x 10 11 cm~ 2 from the O I* A1304.8 A line, 
and to N(Si n*) < 2.6 x 10 12 cur 2 from the 
Si II* A1533.4 A line using the Keck telescope 
spectrum of Q0347-3819 obtained by Prochaska 
& Wolfe (1999). These data are included in Ta- 
ble 3. The most probable b parameters, corre- 
sponding to the column densities from Table 3, 



are as follows: &hi+Mc = 2.71 ± 0.05 km s _1 , and 
bcrii+Fdi = 2.80 ± 0.26 km s _1 . It shows that 
in spite of the difference in shapes of p(v) for two 
groups of ions, the broadening parameter in the 
central H2-bearing component is the same. 

The measured value of log A(H I) = 20.626 ± 
0.005 cm -2 is in good agreement with the DDM 
measurement of log TV (H I) = 20.63 ± 0.09 cm" 2 . 
The re-calculated value of D/H = (3.75 ± 0.25) x 
10~ 5 is consistent with standard big bang nucle- 
osynthesis predictions (e.g. Buries et al. 1999) 
if the baryon-to-photon ratio, rj, is in the range 
4.37 x 10~ 10 < i] < 5.32 x 10" 10 . With the 
present-day photon density determined from the 
CMBR (Mather et al. 1999), we can estimate 
the present-day baryon density fib h1 00 = 3.73 x 
10 7 (T C mbr/2.75) 3 77 ~ 0.016 + 0.020. 

Our new D/H ratio is shifted toward the higher 
hydrogen isotopic ratio found in the Lyman Limit 
Systems (LLS) by Tytler et al. (2000), D/H = 
(3.5 ± 0.5) x 10~ 5 . We recall that the ratio of 
D/H ~ 4 x 10~ 5 was found from the Tytler et 
al. LLSs using a different mesoturbulcnt method 
by Levshakov, Kegel, & Takahara (1998) and by 
Levshakov, Tytler, & Buries (2000). We thus sum- 
marize that when all most accurate measurements 
of D/H in high redshift QSO absorbers are con- 
sidered, one finds the same D/H value which is 
about 2.5 times the mean ISM value of D/H = 
(1.50 ±0.10) x 10- 5 (e.g. Linsky 1998). 

To conclude this section, we comment on the 
fitting of the C II A1036 A line: the exact mea- 
surement of the column density for this ion is very 
important since the calculation of the volumetric 
gas density Tin in the z a b s = 3.025 molecular cloud 
is based on the ratio of N(C II*)/A(C II). The 
C II A1036 A line is saturated and having this 
line alone we would be able to estimate only a 
lower boundary for N(C II). But the system un- 
der study exhibits a plenty of other metal profiles 
which are unsaturated. Being fitted together with 
C II, these profiles allow us to restrict the shape 
of the radial velocity distribution p(v) and then 
to determine the real column density of C II. As 
already stated above, the measurement of N(C II) 
from the saturated A1036 A line is possible if the 
density weighted velocity distribution functions of 
C II and a match ion (which shows unsaturated 
absorption profiles) are equal. 

Indeed, for a fixed value of A (or v in the veloc- 
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i(A) = k N ail / n H (a;) 



ity space) within the line profile, the optical depth 
T a ,i(A) of element 'a' in the ith ionization stage is 
given by [cf. eqs. (24)-(26) in Levshakov, Aga- 
fonova, & Kegel 2000] : 

T a,i[Ef(3Q] 

T 

$ A [AA D (x),«(x)]dx , (15) 

where fc is a constant for a particular line, N a ^ the 
column density, nu{x) the volumetric gas density 
at point x, Y aj , = n a ,i/n a the fractional ionization, 
f = « P h/«H the ionization parameter with n p h 
being the number of photons with energies above 
one Rydberg, <f>\ the profile function at the local 
velocity v(x), and T ayi the mean density- weighted 
fractional ionization. 

If the line broadening is caused by bulk motions 
and the ratio of the fractional ionizations for a pair 
of ions 'a,i', 'b,j' remains the same for all U(x) 
along the sightline, then 



Ta,i(A) N, 



oc 



T b,jW N K 



constant . 



(16) 



Given the integral part of (15), which can be es- 
timated from the unsaturated lines, we can easily 
calculate the column density for even saturated 
lines. 

To test how accurately C II traces the distribu- 
tions of other elements, we have studied the ion- 
ization properties of the warm H I gas assuming 
that it is embedded in an ionizing radiation field 
of given spectrum and intensity. 

The CLOUDY code (Ferland 1996) was used 
to compute the fractional ionization of ions T a> j 
as a function of the ionization parameter U. For 
a given value of n p h, the fluctuation of U along 
the line of sight represents the reciprocal hydro- 
gen density. We consider two limiting types of ra- 
diation fields: a hard, QSO-dominated spectrum 
representative of the radiation field external to 
the H 2 -bearing cloud at z = 3 (Haardt & Madau 
1996), and a soft, starlight background given by 
Mathis, Mezger, & Panagia (1983). For the cho- 
sen ionizing spectra, n p h = 2.2 x 10~ 5 cm~ 3 and 
2.0 x 10~ 6 cm~ 3 , respectively. Calculations were 
performed for metallicity Z = O.IZq. 

We found that most concordant with C II is 
the ion Si II which traces C II over a rather 
wide range of U. This is illustrated in Figure 15 



where the ratios Tsni/Tcn and Tpeii/Tcn (to 
stress the difference) are shown as functions of 
U . The displayed boundaries of U correspond to 
tih = 0.22 -T- 2.2 crn^ 3 (hard ionizing spectrum) 
and to riH = 0.02 0.2 cm~ 3 (soft ionizing spec- 
trum). Figure 15 shows that Tgiii/Tcn is practi- 
cally constant for both ionizing fields. The behav- 
ior of Tpoii/Tcn is, however, different: in both 
cases we may expect to observe non-similar pro- 
files of the C II and Fe II lines. 

Thus we can conclude from these calcula- 
tion that independently on the background ion- 
izing field the distribution of C II does trace 
the distribution of Si II in diffuse clouds (0.1 
cm~ 3 < 7Jh < 10 cm~ 3 ) and therefore we can use 
the unsaturated Si II profiles to measure A^(C II) 
with a sufficiently high accuracy. In the present 
case we were able to describe with the same veloc- 
ity distribution not only the C II and Si II lines, 
but profiles of the other elements. This fact seems 
to justify the suggestion that the density fluctu- 
ations in the H2 cloud are rather small and thus 
the derived N(C II) value can be considered as 
reliable. 

4. Discussion 

We turn now to the question of the physical 
conditions in the z a b s =3.025 DLA. The measured 
column densities of metals and H 2 can be used to 
estimate the volumetric gas density and the in- 
tensity of the UV radiation field. We begin with 
the metal absorptions. Then we consider if the 
observed H2 distribution over the excited rota- 
tional states can be described with metals in a 
self-consistent way. After that we compare the 
revealed [a-element /iron-peak] ratio and abun- 
dances with those known from different DLAs. 

4.1. Metal absorption from the ground 
and excited states 

The absorption lines from the fine-structure lev- 
els of atoms and ions provide comprehensive infor- 
mation about the physical properties in the DLAs 
(e.g., Bahcall & Wolf 1968; Silva & Viegas 2001). 

In equilibrium, the population ratio of the up- 
per level rii to the lower level ni, in ions with a 
doublet fine structure in the ground state, is given 
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by 
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m 



Qifi 



i2.1 



(17) 



where wi,2 and Qi,2 are the photon absorption 
and the collisional excitation rates, respectively, 
and ^2,1 is the radiative transition probability. 

In the present analysis, we are interested mainly 
in the physical conditions in the gas where both 
the H 2 and the metal absorptions have been de- 
tected. We assume that this absorbing region 
is homogeneous and therefore the number den- 
sity ratios are equal to the corresponding ratios 
of column densities, n 2 jn\ — N 2 /Ni. We neglect 
fluorescence since this indirect excitation process 
has usually lower rates as compared to the direct 
collisional excitations and/or the IR radiative ab- 
sorption rates. The radiative decay and the de- 
excitation rates are also ignored since in our case 
they have low efficiency 

In the main component of the z a bs = 3.025 sys- 
tem the following quantities were estimated from 
the data listed in Table 3: A(H I) = (2.521 ± 
0.035) x 10 20 cm" 2 , N(C II) = (5.05 ±0.28) x 10 15 
cm" 2 , C II*/C II = (4.5 ± 0.3) x 10~ 3 , C I/C II 

< 8 x 10~ 5 , O I*/0 I< 3 x 10~ 5 , and Si H*/Si II 

< 2 x ID" 3 . 

The ground state of the C + ion consists of two 
levels 2s 2 2p 2 ~P\/ 2 3/2 w ^h the energy difference of 
AE = 91.32 K and the spontaneous decay rate of 
A 2 .i = 2.291 x 10~ 6 s" 1 . The excited J = 3/2 level 
may be populated by collisions with electrons, e~, 
hydrogen atoms, H°, protons, H + , and molecules, 
H2 . Among them the highest rate is found for col- 
lisions with electrons : Q\ 2 — 2 x 10~ 7 n c s -1 in 
the range 10 2 K < T kin < 10 3 K (see, e.g., Figure 3 
in Silva & Viegas 2001). But the electron density 
in H I regions comes mainly from the photoioniza- 
tion of carbon, n c ~ (C/H)tih, which is approx- 
imately equal to 3 x 10~ 5 rin for the z a b s = 3.025 
DLA. Then Q\ 2 ~ 6 x 10" 12 n H s" 1 becomes 
lower than the H°-C + collisional rate for the same 
temperature interval where Q^ 2 is almost constant 
and equals Qf° 2 — 1-45 x 10~ 9 «h s _1 (Launay & 
Roueff 1977a). The rate of the photon absorption 
in the z a b s = 3.025 DLA is mainly induced by the 
cosmic microwave background radiation (Molaro 
ct al. 2001b) and equals to wi, 2 — 1.1 x 10~ 9 s _1 . 

From (17) and taking the limiting values of the 
measured C II* /C II ratio, one obtains rh — 6 -f- 7 



cm 3 . The low limit of tir 4 cm 3 stems from 
the analysis of the H 2 lines (see §4.2). 

To check the consistency of the estimated tin 
with other metal absorption lines, we can use the 
upper limits to the column densities of O I* and 
Si II* . The equilibrium equation for the three- level 
system ( 3 P2, 3 Pi, 3 Po states of O ) can be written 
as (cf. Molaro et al. 2001a) 



m (<$°i + Q£° 2 ) 

(NJ I - 

n ~ Mfi 



(18) 



where the subscripts 0, 1, and 2 denote the ground 
(0) and excited (1,2) levels, respectively. The 
spontaneous decay rate Ai, = 8.95 x 10~ 5 s _1 , 
and the energy differences between the ground 
state and the corresponding fine-structure levels 
are AE ,i = 228 K and AE . 2 = 326 K. 

If n c is small, the rate coefficient for collisions 
with hydrogen atoms is dominating. At the kinetic 
temperature Tkin = 400 K, we find from Launay 
& Roueff (1977b) Q£i = 1.65 x 10" 11 n H s" 1 and 

Q£° 2 = 0.10xl0- n n H s" 1 . Then (18) yields n H < 
150 cm~ 3 , which is consistent with the foregoing 
value. 

The fine-structure level of the Si + ion { 3 P^i 2 ) 
has an effective excitation temperature ~ 414 K 
which means that the estimated upper limit to 
A (Si II*) is less sensitive to low gas densities in 
this last case. Besides, absorption features pro- 
duced by Si II* must come only from H II regions 
as argued by Spitzer & Jenkins (1975). The ab- 
sence of the Si II* absorption might indicate that 
there is no dense H II gas region on the line of 
sight. This suggestion is also supported by the ab- 
sence of the N II* A1084.562 A and A1084.580 A 
lines in the £ a b s = 3.025 system. Both facts show 
that the C II* and C II absorptions arise predom- 
inantly in a warm H I gas in agreement with our 
model described in §3.6. 

4.2. H 2 distribution over the rotational 

states 

We find no pronounced differences between the 
H 2 abundances in the z a b s = 3.025 DLA system 
and in the Galactic diffuse clouds which have com- 
parable H I column densities. The high excitation 
temperature in our system, T cx = 825 ± 110 K, 
is also consistent with previous observations in 
the Galaxy. For instance, in a pioneer work by 
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Spitzer & Cochran (1973) it has been shown that 
T cx > 300 K within diffuse interstellar clouds 
where the H2 column density in the J = rota- 
tional level iV(0) < 10 15 cm" 2 , and the tempera- 
ture as high as 1100 K has been detected by Mor- 
ton & Dinerstein (1976) toward £ Pup. 

It has been discussed in §3.4 that the popula- 
tions in the excited rotational states of H 2 cannot 
be supported only by collisions in a gas whose ki- 
netic temperature is as high as 800 K. We cal- 
culated the H 2 collisional excitation rates by us- 
ing the analytical approximations given in Shull 
& Beckwith (1982). The deduced rates are listed 
in Table 5, columns 3 and 4 for T kin = 400 K 
and 1000 K, respectively. In columns 5 and 6 we 
show, for comparison, the rates obtained in dif- 
ferent models by Nishimura (1968) and Elitzur & 
Watson (1978). As noted by Shull & Beckwith, 
the accuracy of the collisional excitation rates are 
not sufficiently high and different models may pro- 
duce deviations of about 50%. The employed an- 
alytical fit to the H°-H2 excitation rates agrees 
with the Elitzur & Watson calculations to better 
than 25 %, but gives much lower values than the 
Nishimura rates. The radiative rates listed in col- 
umn 2 were taken from Turner, Kirby-Docken, & 
Dalgarno (1977). 

The collisional excitation of the J = 3 level is 
not expected to be effective in our case. Indeed, 
the J = 2 level has longer radiative lifetime when 
compared with the J — 3 level (tj=2 = 3.4 x 10 10 
s and tj = 3 = 2.1 x 10 9 s, respectively), and there- 
fore the J = 2 level is more sensitive to the col- 
lisional de-excitation. The critical density, = 
^2^0/92^0: above which collisional de-excitation 
becomes important is = 3.6 cm~ 3 . Thus, 
nn > 4 cm~ 3 is required to maintain the observed 
N(2)/N(0) ratio at T kin = 400 K. The absence of a 
fracture at lower J in the graph log [AT (J) /g(J)] ver- 
sus E(3), shown in Figure 8, can be used to infer 
an upper limit for the gas density, 7Jh J$ 74 cm~ 3 . 
Both limits are consistent with «h estimated from 
the measured N(C ll*)/N(C II) ratio. 

Now we can set restrictions on the rate of 
photo-dissociation, /, in the z a b s = 3.025 molec- 
ular cloud assuming that the H 2 molecules are 
formed on grains with rate R and partly destroyed 
through the fluorescence process. In equilibrium, 
the fractional molecular abundance is given by 



(Jura 1974) 

/h 2 ~ 2Rn H /I ■ (19) 

To estimate the rate of formation of H2 , we use 
the approximation from Black & Dalgarno (1977) 

R ~ 10~ 6 T k ^ n H n d cm 3 s^ 1 . (20) 

Here the number density of dust grains, n<j, is re- 
lated to the total hydrogen density, n, through 

n d = 10- 12 n = 10- 12 (n H + 2n H2 ), (21) 

which is the usual relation for the interstellar dif- 
fuse clouds. 

Then, accounting for the dust-to-gas ratio k ~ 
0.05 (see § 3.5) and the fact that riH 2 <C nu in the 
^abs = 3.025 cloud, we find 

R ~ 10- 18 T k ° in 5 k n\ ~ 5 x 10- 20 T k ° in 5 n\ cm 3 s" 1 . 

(22) 

Then for T kin = 400 K, H 2 forms at a rate R ~ 
1.0 x lO" 18 ^ cm 3 s" 1 . 

This value of R will be consistent with the up- 
per limit found by Jura (1974) in the Galactic 
ISM, Rq < 2xl0~ 16 cm 3 s" 1 , ifn H < 14 cm" 3 at 
z a bs = 3.025. By this we confine nn in the interval 
4 cm -3 < nn 5; 14 cm~ 3 which provides a self- 
consistency between the observed N(2) /N(0) ratio 
and the maximum formation rate of H2 under the 
condition that the physical properties of the dust 
grains in the Milky Way and at z a b s = 3.025 are 
similar. 

Since the molecular hydrogen fractional abun- 
dance in the main component of the z a b s = 3.025 
absorber is /h 2 = (3.3 ± 0.2) x 10~ 6 , cloud pa- 
rameters nn — 6 cm" 3 and T k i n = 400 K will 
yield / ~ 2 x 10~ 10 s _1 which agrees well with 
the values measured in the interstellar clouds con- 
taining optically thin H 2 : 0.5 x 10~ 10 s" 1 < I G < 
16.0 x lO" 10 s" 1 (Jura 1975a). Thus, the radiation 
field in the z a b s = 3.025 cloud is quite compara- 
ble with those measured in transparent interstellar 
clouds in the MW. 

We now evaluate the total hydrogen density, 
nH, the photo- absorption rate, f3o, and the steady 
state production rate of molecular hydrogen, i?, 
from the relative population of the different ro- 
tational levels of H 2 . Our analysis is, in gen- 
eral, the same as described by Spitzer & Zweibel 
(1974) and by Jura (1974, 1975a,b). However, 
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we re-calculated the formation efficiency coeffi- 
cients £4 and £5 for the J = 4 and 5 levels in- 
cluding upper rotational levels with a cutoff above 
Jmax = 15. These coefficients were calculated from 
the Boltzmann distribution of molecules formed 
upon grain surfaces with an effective tempera- 
ture, Tf (for details, see Black & Dalgarno 1976). 
Since a newly formed molecule quickly cascades 
to the ground vibrational level, we assume that 
the newly formed molecules populate the (v=0, 
J=4) and (v=0, J=5) levels. It is usually as- 
sumed that molecules are formed at Tf ~ 17400 K, 
which is the temperature equivalent of the bind- 
ing energy of H 2 distributed statistically among 
the rotation-vibration levels upon formation (e.g. 
Spitzer & Zweibcl 1974). Then we find e 4 = 0.21 
and £5 = 0.65 (cf. £4 = 0.19 and £ 5 = 0.44 from 
Jura 1975b). 

The concentrations of H2 in the 3 — 4 and 5 
levels may be approximately estimated from the 
statistical equilibrium equations 

^4+K,2+K «4+K = £4+K Rn\\ + (3q ^ PA+K..T nj> , 

J' 

(23) 

where J' = K, K + 2, K + 4, and K = or 1 for J 
= 4 or 5, respectively. 

In these equations, the pumping efficiency co- 
efficients pi j are taken from Jura (1975a). The 
quantity pij shows the probability (averaged over 
the different Lyman and Werner bands) that 
photo-absorption from the jth rotational level will 
lead to the population of the ith rotational level 
in the ground v = vibrational state. 

From (23) we can calculate Run using the mea- 
sured iV(4)/iV(H) and N(5)/N(R) ratios. It gives 
Rn H ~ 3 x 10~ 16 s" 1 . If we take R ~ 5 x 10~ 17 
cm 3 s _1 from Spitzer & Cochran (1973), then 
riH — 6 cm~ 3 . Now the H 2 destruction rate 
can be estimated from (19) which gives again 
/ ~ 2 x 10~ 10 s _1 as in the foregoing estima- 
tions. If I ~ 0.11/? (Jura 1974), then the rate 
of photo-absorption is /3o ~ 2 x 10~ 9 s _1 . On the 
other hand, the average flo for the intergalactic UV 
background at z ~ 3 is about 2 x 10~ 12 s _1 (Sri- 
anand & Pctitjean 1998), while the MW interstel- 
lar radiation field value ranges between 5 x 10~ 10 
< A) < 16 x 10~ 9 s- 1 (Jura 1975a). Thus, our 
measurements show a close similarity with the av- 
erage characteristics of the UV field found in the 



MW. 

For the neutral hydrogen column density 
N(R I) ~ 2.5 x 10 20 cur 2 and n H ^ 6 cm" 3 , 
the linear size of the H2 -bearing cloud along the 
line of sight is equal to D ~ 14 pc which yields 
a mass of M ~ 0.5D 3 nnmn ~ 2OOM , assuming 
spherical geometry. 

4.3. [a-element /iron-peak] ratios 

The elemental column densities obtained fol- 
lowing the procedure described in sections 3.5 and 
3.6 together with the relative abundances are sum- 
marized in Table 4. 

Among the elements belonging to the iron-peak 
group, Zn is a good indicator of the absolute metal 
abundance in the system, since Zn is known to 
be essentially undepleted in the ISM (Pettini et 
al. 1999). Measurements of Zn in high rcdshift 
systems have been proven to be very challenging 
since the Zn II resonance lines fall at wavelengths 
> 8000 A where spectrograph and CCD efficiency 
decrease while sky emissivity increases. The Zn 
measurement presented here and that in the DLA 
at z a bs = 3.39 toward Q0000-2621 are the only 
ones obtained for systems with z > 3. The derived 
[Zn/H] = -1.50 ± 0.05 is in line with the bulk of 
[Zn/H] measurements found in the redshift inter- 
val 2 < z < 2.5 although higher by 0.6 dex than 
that of the DLA at z abs =3.39 toward Q0000- 
2621 confirming the presence of significant spread 
in absolute abundance in DLA with comparable 
redshifts. The straight mean of the Zn measure- 
ments in the redshift interval from 2 < z < 3 is 
[Zn/H] w —1.2, and does not differ from the col- 
umn density weighted mean for the same redshift 
interval (Vladilo et al. 2000). This new measure- 
ment therefore supports the general metallicity be- 
havior, characterized by the intrinsic spread for a 
given rcdshift and by a mild trend of increasing 
abundance with the decreasing redshift if not a 
lack of evolution when the hydrogen column den- 
sity weighted mean is considered. 

As described in previous § 3.6, our model keeps 
the constant ratio between the matching compo- 
nents. With this condition imposed, we cannot 
account for different chemical abundances or the 
presence of different dust pattern depletion among 
the components of the system. It may be possible 
that a contribution from warmer gas or different 
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dust depletion is responsible for the observed dif- 
ference in profile shapes of Fe II and other ions. 

Only recently it became possible to measure 
oxygen with an acceptable accuracy using the 
fainter O I lines falling in the Lya forest instead of 
the always strongly saturated O IA1302 A tran- 
sition (Lopez et al. 1999; Molaro et al. 2000; 
Dessauges-Zavadsky et al. 2001). In the present 
paper, the O abundance is derived from the set of 
eleven OI lines used either for direct fitting or just 
for the consistency test. 

The [O/Zn] is the most informative ratio to 
trace back the kind of chemical evolution of the 
absorber. Both elements show little affinity with 
dust grains and are good representatives of en- 
chain and iron-peak elements, respectively. If the 
system is in the early stages of chemical evolution, 
the relative abundances should reflect the yields of 
Type II SNe characterized by the enhancement of 
en-elements. In fact, it was rather surprising to 
find little evidence for this enhancement when el- 
ements such as S and Si, once corrected by dust, 
where used in combination with Zn (Centurion et 
al. 2000; Vladilo 1998) and the ratio [O/Zn] = 0.1 
has been measured in the DLA toward Q0000- 
2621 (Molaro et al. 2001a). 

The ratio [O/Zn] = 0.68 ± 0.08, found here, 
represents the first clear cut evidence for a gen- 
uine enhancement of the a-elements comparable 
to that observed in the metal poor stellar popu- 
lation of the Milky Way (Israelian et al. 2001). 
The oxygen result is also supported by the other 
a-elements we have been able to measure, namely 
[Ar/Zn] = 0.34 ± 0.07, [Si/Zn] = 0.56 ± 0.06, and 
[S/Zn] = 0.41 ± 0.08 (if we take logiV(S II) = 
14.73 ± 0.01 cm" 2 from Prochaska & Wolfe 1999). 
Argon abundance is slightly below the values ob- 
tained from the other a-elements, but this is not 
surprising since Ar is sensitive to photoionization 
effects and shows in the Galaxy even larger deple- 
tion factors (Sofia & Jenkins 1998). 

All results available at present allow us to con- 
clude that the DLAs probably comprise a mix of 
chemical abundance pattern with both low and 
standard a-element enhancement. It is not clear 
at the moment if this reflects the different enrich- 
ment histories or just the same type of galaxies 
caught at different stages of their chemical evolu- 
tion. It is interesting to note that estimated in 
the H 2 -bearing cloud the ratio (Ar/O) = —2.65 



is comparable to that observed in the DLA to- 
ward Q0000-2621, (Ar/O) = -2.25 (Molaro et al. 
2001a), in particular when allowing for a small de- 
gree of Ar photoionization. As argued by Henry & 
Worthcy (1999) the similarity of these ratios in a 
variety of galactic systems from the MW to Blue 
Compact galaxies strengthens the case for the ex- 
istence of a universal IFM. 

It is commonly assumed that nitrogen is 
mostly produced by the intermediate massive stars 
(Henry, Edmunds, & Koppen 2000). In DLAs, N 
shows a scatter well in excess of observational un- 
certainties which has been interpreted as being 
due to the delayed release of N compared to O 
(Centurion et al. 1998). Our value in the '[O/H]- 
[N/O]' plane has coordinates [O/H] = -0.82 and 
[N/O] = —0.84 and falls close to the line for a pure 
secondary origin for N (cf. Figure 8 in Centurion 
et al. 1998). If intermediate massive stars had 
not yet time to evolve and to contribute to the 
N enrichment then this sets an upper limit of ap- 
proximately 10 8 yr to the evolutionary age of the 
z a bs = 3.025 system. The young age of the system 
under study is also consistent with the observed 
en-element enhancement. 

Phosphorus has been measured in the DLA to- 
ward Q0000-2621 (Molaro et al. 2001a) and to- 
ward Q1759+7539 (Outram et al. 2000). P cannot 
be measured in halo stars and thus the damped 
Lya systems provide a unique astrophysical site 
where this element can be measured in low metal- 
licity environments. Our determination is based 
only on the P II A963.8 A line and may be slightly 
overestimated due to H I contamination. We have 
here the mean value of [P/Si] « —0.2 which sug- 
gests the presence of a moderately odd-Z/even-Z 
effect which is a common feature in all P measure- 
ments in DLAs so far. 

5. Conclusions 

The main results of this paper are as follows. 

(1) In the high resolution UVES spectrum of 
the quasar 0347-3819 we have identified over 80 
absorption features with the Lyman and Werner 
lines arising from the J = to 5 rotational levels 
of the ground electronic-vibrational state of H 2 . 
These lines belong to the z a b s = 3.025 damped Ly- 
a system. 

(2) The application of the standard Voigt fit- 
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ting analysis to the main subcomponent at zjj 2 = 
3.024855±0.000005 gives the total H 2 column den- 
sity of iV(H 2 ) = (4.10 ±0.21) x 10 14 cm" 2 and the 
Doppler width &h 2 = 2.80 ± 0.45 km s _1 . The 
kinetic gas temperature is less than 430 K. 

(3) The fractional abundance of H2 in the main 
component is equal to /h 2 = (3.3 ± 0.2) x 10~ 6 , 
while being compared to the total H I column den- 
sity it yields /h 2 = (1.94±0.10) x 10~ 6 . These val- 
ues are similar to those observed in the ISM diffuse 
clouds with low color excesses, E(B — V) < 0.1. 
Using the total hydrogen column density in the 
main component N(H I) ~ 2.5 x 10 20 cm~ 2 and the 
approximate relationship N(R) ~ 7.5 x 10 21 E(B — 
V) cm~ 2 mag -1 , which assumes that the gas and 
dust grains are uniformly mixed (Spitzer & Jenk- 
ins 1975), we obtain E(B - V) ~ 0.03 mag in the 
H 2 -bearing cloud. For the second strong compo- 
nent with JV(H I) ~ 1.5 x 10 20 cm -2 at Av ~ -17 
km s _1 we set a very low limit to /h 2 < 3.3 x 10 -9 
(3cr). The z a bs~ 3.025 system is not a homoge- 
neous cloud but rather a mixture of cold and warm 
gas. 

(4) The total column densities in the ortho- 
and para-H 2 states yield the ratio -/V ort ho/^para = 
3.1 ± 0.3 which is close to 3:1 expected from the 
hot H 2 formation upon grain surfaces. Similar 
ratios of 2.4 ± 0.9 and 3.1 ± 0.8 are estimated 
at z abs = 2.33771 and 2.8112 toward Q1232+0815 
and Q0528-250 from the published data by Sri- 
anand et al. (2001) and Srianand & Petitjean 
(1998), respectively. These ratios are, however, 
higher than the 'freeze-out' ortho:para-H 2 =1:4 
predicted for the early universe at z < 20 (Flower 
et al. 2000). It appears that H 2 detected in these 
three DLAs has undergone a re-forming phase in 
the period between z ~ 20 and z ~ 3 (At ~ 2 
Gyr). 

(5) We found that the population of the low 
rotational levels can be represented by a single ex- 
citation temperature of T ex = 825 ± 110 K. Such 
high T ex agrees well with the results by Spitzer 
& Cochran (1973) who showed that T cx > 300 K 
within diffuse clouds which have N(3 = 0) < 10 15 
cm~ 2 , as it is in our case. 

(6) From the measured [Cr/H] and [Zn/H] ra- 
tios we calculated the gas-to-dust ratio k = 0.032± 
0.005 in the whole z a b s = 3.025 system (the ra- 
tio [Cr/Zn] is usually used to indicate the pres- 
ence of dust in DLAs). The derived H 2 fractional 



abundance fits well in the linear relation between 
log/n 2 and [Cr/Zn] supporting our previously ob- 
tained results for other DLA systems (Levshakov 
et al. 2000b). 

(7) From the observed population ratios of 
H 2 in different rotational states and the C II* /C II 
ratio we infered the hydrogen number density tih 
in the H 2 -bearing cloud «h — 6 cm~ 3 , and the 
cloud dimension along the line of sight D ~ 14 pc 
which gives M ~ 200M Q . 

(8) From the relative populations of H 2 in the 
J = 4 and 5 rotational levels we estimated the 
rate of photo-destruction / ~2x 10~ 10 s _1 . The 
photo-absorption rate, (3q ~ 7/0.11, is thus equal 
to /3 — 2 x 10~ 9 s _1 . Taking into account that 
the MW interstellar radiation field value ranges in 
5 x 10~ 10 s- 1 < f3 < 16 x 10~ 9 s- 1 (Jura 1975a), 
we conclude that the UV radiation fields in the 
z&bs = 3.025 absorbing cloud and in the Galactic 
ISM are very much alike. We also found that the 
formation rate of H 2 upon grain surfaces, Run ~ 
3xl0~ 16 s _1 , is similar to that found in the Galaxy 
by Jura. 

(9) The metal abundances from the z a b s = 3.025 
DLA reveal a pronounced [a-element/iron-peak] 
enhancement with [0,Si/Zn] = 0.6 ±0.1 at the 6<r 
confidence level, the first time that this abundance 
pattern is unambiguously identified in a DLA sys- 
tem. The same order of magnitude enrichment 
of a-chain elements has been observed in Galac- 
tic metal-poor stars. The measured [N/O] ratio 
implies the chemical history age of this DLA of 
< 10 8 yr. 

(10) The analysis of the metal profiles combined 
with the H+D Lyman series lines yields a new es- 
timation of the hydrogen isotopic ratio D/H = 
(3.75 ± 0.25) x 10~ 5 in this DLA. It implies for 
the baryon-to-photon ratio, 77, a value in the in- 
terval 4.37 x 10- 10 < 77 < 5.32 x lO" 10 , and the 
present-day baryon density fib h 2 ~ 0.016 0.020, 
which is in good agreement with a new anal- 
ysis of the BOOMERANG experiment yielding 
n h h 2 = 0.021 ± 0.003 (Nctterfield et al. 2001). 
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A. Correlated measurements 

In order to evaluate the uncertainties in the derived physical parameters, one must know the standard 
errors, aj, for all the data points. These errors, computed through the UVES data reduction pipeline 
(Ballcster et al. 2000) are not, however, statistically independent. The pipeline data resampling to constant 
wavelength bins introduces a correlation between the data point values. Besides the resampling procedure 
can change statistical properties of the noise (e.g., the simplest linear interpolation scheme smooths the 
random noise since the uncertainty of the data point value half way between two measurements is less than 
that for the real detector amplitudes at the limiting points). 

Data point correlation is clearly seen in Figure 7 where the sequences of points at the continuum level 
(e.g., panels L2-0R0 or L2-0P3) show that neighboring points differ from each other by less than dj (for the 
data shown in Fugure 7, the correlation coefficient r ~ +0.77). The correlation length, £, of about 1.5 bin 
can be deduced from the analysis of the autocorrelation function (ACF) 11 of the intensity fluctuations in 
the continuum windows free from absorption-line features. This ACF is shown by histogram in Figure 16a 
where an exponential correlation function /(An) = cxp(— An/£) of the simplest data model, the first order 
autoregression (the Markov process), - is overdrawn by the smooth curve (here An is the scale shift). 

The correlated data make the standard A\ 2 calculation of parameter confidence intervals questionable. 
To investigate this point we performed the following Monte Carlo (MC) analysis in combination with the 
least-squares method. 

It is well known that the least-squares method (LSM) is computationally equivalent to the ^-minimization 
(e.g., Meyer 1975). It is also known that the LSM is an unbiased and minimum variance estimator (Gauss- 
Markov theorem) irrespective of the distribution of the measurements including correlated data as well 
(Aitken 1934). Therefore, in the MC analysis, we used the best fitting parameters (the parameter vector 9, 
see § 3.6) found by means of the standard x 2 -minimization, and calculated uncertainties of the ^-components 
in the vicinity of the global minimum of the objective function by using the MC procedure. 

Our computational MC procedure consists of the following steps: 

(1) A simulation box in the parameter space is specified by fixing the parameter boundaries 9i(l — (3) < 
9i < 6i(l + /?); in our case f3 = 0.3. 

(2) Initial guess 9^ is chosen arbitrary within the simulation box. 

(3) A random realization of the 'observational' points is generated in the followin way: the jth data point 
value is altered with a probability P by means of the relation 

J-j = J- 3 + <Jj g , 

where g is a random normal variable with zero mean and unit dispersion, and J-° hs is the normalized 
observational intensity 12 in the jth bin (j = 1, 2, . . . , A/"; with M being the total number of data points). 

(4) Objective function is minimized and all solutions with x 2 < 1-10 are collected. 

(5) The whole procedure is repeated many times to collect a representative sample (i.e. the number of 
measurements K. > 30). Then the sample mean and standard deviation are calculated. 

To compare both A% 2 and MC results we carried out the MC procedure for the b parameter and the 
H 2 abundances at different J levels. Nineteen H 2 lines (L3-0R0, L7-0R0, L1-0P1, L3-0P1, L14-0P1, W0-0R2, 
W0-0Q2, L3-0R2, L4-0R2, L5-0R2, L2-0R3, L2-0P3, L3-0P3, L4-0R3, L15-0R3, L5-0R4, L8-0P4, L7-0P5, 
L12-0R5; see Figs. 5 and 6) have been chosen and fitted simultaneously with the individual Gaussian profiles. 
The obtained estimations (for M = 177 and K. = 181) summarized in Table 6 show that in this case Ax 2 yields 

11 The ACF is defined, e.g., in Bendat & Piersol (1971) 

12 This randomization decreases the correlation between consecutive bins, but increases the noise. For instance, if P = 0.7 the 
correlation length equals £ ~ 1.0 bin (see Figure 16b) but the noise level is about a factor of \/2 higher. 



19 



practically the same mean values but systematically higher standard deviations. Since statistical properties 
of the noise are not exactly known we use in the present paper the most conservative A% 2 estimations of the 
standard deviations and of the mean values. 



B. Metal systems toward QSO 0347 3819 

We have thoroughly investigated all possible coincidences of H 2 , H I, D I, and metal absorption features 
from the z a b s = 3.025 damped Ly-a system with metal absorptions from other systems observed toward 
Q0347-3819. Positions of the found blends are indicated in Figures 5, 6, 10, and 11. The strongest systems 
detected in the spectral ranges AA = 3680 — 4880 A. 4880 — G730 A (this portion of the spectrum obatined 
at the Keck telescope was kindly shared with us by Wolfe & Prochaska), 6730 - 8520 A, and 8685 - 10000 A 
are listed below. Since almost all identified metal profiles exhibit complex structures, the quoted redshifts 
are the ones of the strongest components or the average redshifts in the cases when the components are not 
well separated. In cases where the hydrogen Ly-a line lies outside the observable spectral range, we consider 
a doublet of absorption lines caused by the same ion as an individual metal system if the lines show identical 
profile shapes. The identifications made in the z a b s = 3.025 DLA are not included in the list. Possible blends 
are shown in parenthesis. 



21 = 


= 0.000019 


Ca II 3934,3969(+Wl-0R2) 


22 = 


= 1.4579 


Mg I 2852; Mg II 2796,2803; Al II 1670; Al III 1854; Si II 1526,1808; 






Fe II 1608,2260,2344,2374,2382,2586,2600 


23 = 


= 1.5159 


CIV 1548,1550 


24 = 


= 1.5263 


Mgll 2796,2803; Sill 1526(+Ly- 7 2i 5 ),1808; Fell 2344,2382,2600 


25 = 


= 1.6768 


CIV 1548,1550 


26 = 


= 2.2140 


Ly-a; C II 1334; C IV 1548,1550; O I 1302; Mg II 2796,2803; Al II 1670; 






A1III 1854; Sill 1190(+W3-0Q3), 1193,1260; Si III 1206; Si IV 1393,1402; 






Fe II 1608,2344,2374,2382,2586,2600 


27 = 


= 2.3932 


Ly-a(+Ly-/3 DLA); C IV 1548,1550 


28 = 


= 2.4344 


Ly-a; C IV 1548,1550 


29 = 


= 2.5371 


Ly-a; C II 1334; C IV 1548,1550; Al II 1670; Si II 1260; Si IV 1393 


210 


= 2.6504 


Ly-a; C IV 1548,1550; Si III 1206; Si IV 1393,1402 


211 


= 2.6526 


Ly-a; C IV 1548,1550; Si III 1206; Si IV 1393,1402 


212 


= 2.8102 


Ly-a; C II 1334; Sill 1190, 1193, 1260(+Si II 1193 DLA); Si IV 1393,1402 


213 


= 2.8483 


Ly-a, Ly-/3(+Ly-7 2i 9 ); C III 977; C IV 1548,1550 


214 


= 2.8992 


Ly-a, Ly-/3, Ly-7, Ly-8; C IV 1548,1550; C III 977; 






Si III 1206; Si IV 1393,1402 


215 


= 2.9617 


Ly-a, Ly-/3, Ly-7, Ly-8; C III 977; C IV 1548,1550(+Si II 1526 DLA); 






Nil 1083; NIII 989; Al III 1854; Sill 1190; Si III 1206 


216 


= 2.9659 


Ly-a, Ly-/3, Ly- 7 (+Si II 1526 24), Ly-8, Ly-e; C III 977; C IV 1548; 






NIII 989(+Lll-0P2); Si III 1206 


217 


= 2.9792 


Ly-a, Ly-/3, Ly-7; C II 1036(+Ly-/3 DLA); C III 977; C IV 1548,1550; 


218 


= 3.0217 


Ly-a(+Ly-a DLA); C III 977(+0 I 976 DLA); C IV 1548,1550; N II 1083; 






N III 989; Si IV 1393,1402 


219 


= 3.0632 


Ly-a, Ly-/3, Ly--y (+Ly- f3 z 13 ), Ly-<5(+N III 989 214); C III 977; C IV 1548,1550 
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Fig. 1. — A portion of the normalized spectrum of 
Q0347-3819 in the region of a few H2 lines from 
the Lyman and Werner bands (labels 'L' and 'W, 
respectively) identified in the z a b s = 3.025 damped 
Lya system. The noise level is indicated by the 
dashed curve at the bottom of the panel. The 
general continuum was estimated by using a poly- 
nomial fit to the 'continuum windows' in all spec- 
tral range spanning about 1200 A. In cases where 
the curvature of the QSO continuum at z — 3.025 
is strong, the uncertainty in the local continuum 
level may be dominated by uncertainties in the 
fitting procedure. 

Fig. 2. — Confidence regions in the 'x 2 — V 
plane calculated from the simultaneous fit of the 
H 2 lines from Table 1 as described in the text. The 
parabola vertex corresponds to the most probable 
value of b = 2.80 km s" 1 . 

Fig. 3. — Curves of growth for 31 H2 lines (J = 
0, 1, 2, 3, 4, and 5; see Table 1) in the z a b s = 
3.025 molecular cloud toward Q0347-3819. Most 
probable Doppler parameter is 6jj 2 = 3 km s . 
Curves with 6 = 2 km s _1 and 4 km s _1 restrict 
the lcr confidence region for &h 2 • 

Fig. 4. — Confidence regions in the 'x 2 — LogiV(J)' 
planes calculated from the simultaneous fits of 
the H 2 lines arising from the same J levels. The 
parabola vertexes correspond to the most proba- 
ble values of JV(0) = 1.90 x 10 13 cm" 2 , iV(l) = 
1.60 x 10 14 cm" 2 , N(2) = 5.40 x 10 13 cur 2 , 
N(3) = 1.22 x 10 14 cm" 2 , JV(4) = 2.60 x 10 13 
cm" 2 , and AT(5) = 2.90 x 10 13 cur 2 . 

Fig. 5. — H 2 absorption features associated with 
the z a b s = 3.025 damped Lya system toward 
Q0347-3819 (normalized intensities are shown by 
dots with 1(7 error bars). The zero radial velocity 
is fixed at zh 2 = 3.024855. Smooth lines are the 
synthetic H 2 spectra calculated for the mean phys- 
ical parameters (see text). Possible blends with 
lines from the z a b s = 3.025 system are shown with- 
out indicating redshift, whereas blends with lines 
from other systems are indicated by ionic species 
with redshifts in accord with the list in Appendix 
B. H 2 profiles disturbed by the Ly-a forest absorp- 
tion are also included. 



Fig. 6. — Same as Figure 5. 

Fig. 7. — Spectral regions including the 
H 2 transitions of the Lyman and Werner bands 
(labeled in panels) in the z — 3.024637 absorber 
in the direction of Q0347-3819. The lowest panel 
presents the result of a stacking. The 3<7 upper 
limit of 0.6 mA was calculated for a 12 km s _1 
window centered at v — km s _1 . 

Fig. 8. — Population of the low rotational lev- 
els of H 2 in the z a b s = 3.025 damped Lya sys- 
tem toward Q0347-3819. Error bars for lcr de- 
viations are shown. The negative inverse slope of 
the best-fit straight dotted line indicates the exci- 
tation temperature of about 800 K. 

Fig. 9. -— Relation between H 2 fractional abun- 
dance /h 2 plotted on a logarithmic scale and rel- 
ative heavy element depletion in DLAs. The solid 
line corresponds to the linear law: log/n 2 — 
—6.1 [Cr/Zn] — 7.9, obtained for the five measure- 
ments of /h 2 we know (DLAs # 5, 8, 12, 13 and 
14). The filled circles correspond to [Cr/Zn] and 
the open ones to [Fe/Zn]. See Table 2 for refer- 
ences. 

Fig. 10. — Metal absorption features associated 
with the z abs = 3.025 DLA toward Q0347-3819 
(normalized intensities are shown by dots with la- 
error bars). The zero radial velocity is fixed at 
Zabs = 3.024855. The smooth lines are the syn- 
thetic spectra convolved with the instrumental res- 
olution of FWHM = 7.0 km s _1 . Two groups 
of lines (H I, C II, C II*, N I, O I, Si II, P II, 
Ar I, Zn II, and Fe II, Cr II) are fitted simulta- 
neously to the observed profiles or to their por- 
tions marked by horizontal lines which also indi- 
cate pixels involved in the optimization procedure 
(the corresponding hydrogen profiles are shown in 
Figures 11 and 12). For Zn II and Cr II composite 
stacked spectra are shown. Two absorption fea- 
tures at Av ~ 17 km s -1 (Zn II) and Aw ~ 6 
km s _1 (Cr II) are telluric lines. Possible blends 
with lines from other systems are indicated by 
ionic species with redshifts in accord with the list 
in Appendix B. 
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Fig. 11. — Same as Figure 10 but for the hydrogen 
and deuterium lines. For the Ly-a line only a por- 
tion of its blue wing free of the forest absorptions is 
shown. The blue wing of Ly-/3 is contaminated by 
the forest absorption. The deuterium absorption 
seen in Ly-8, Ly-10 and Ly-12 is shown zoomed in 
Figure 12. The corresponding D/H ratio is equal 
to 3.75 x 1CT 5 . 

Fig. 12. — Same as Figure 11 to show details of 
the H+D (thick lines) and D (thin lines) profiles. 

Fig. 13. — Radial velocity distributions p(v) for 
the models adopted in the present paper (solid 
and dotted curves) and, for comparison, p(v) from 
DDM (dashed curve). The difference in shapes 
of p(v) leads to the difference in the estima- 
tions of D/H: (2.24 ± 0.67) x 10~ 5 (DDM) and 
(3.75 ± 0.25) x 10" 5 (see Figure 14). The intense 
and narrow spike at Aw ~ km s _1 was revealed 
from the analysis of the H 2 lines. 

Fig. 14. — Confidence regions in the 'HI - D/H' 
plane calculated from the simultaneous fit of the 
H I and metal lines shown in Figures 10 and 
11 when the other best-fitting parameters of the 
adopted model are fixed. The contours represent 
68.3% (inner) and 95.4% (outer) confidence levels. 
The cross marks the point of maximum likelihood 
for the adopted model (x^in = l- 22 > v = 545). 

Fig. 15. — Fractional ionization ratios Tsni/Tcn 
and Tpeii/Tcn versus the ionization parameter 
U in the case of a soft, starlight background (solid 
curves) and a hard, QSO-dominated ionizing field 
(dashed curves). Similarity of the line profiles of a 
pair of ions requires the fractional ionization ratio 
to be constant for different U . 

Fig. 16. — Autocorrelation functions (ACF) of the 
intensity fluctuations in the continuum between 
A 6790 A and A 6820 A in the UVES spectrum of 
Q0347-3819 before and after randomization (his- 
tograms in panels a and b, respectively). The 
smooth curves show the ACFs of the first order 
autoregression process with correlation lengths in- 
dicated by £ (see text for more details). 



Table 1 

H 2 Lines at z abs = 3.025 toward Q0347-3819 



J 


Line 


(A), A 


CT(A), A 


W\, mA 


ow x , mA 




z 




t7 z 





L3-0R" 


4277.96 


0. 


09 


6.7 


4.4 


3. 


02486 


0. 


.00009 




L7-0R a,b 


4076.44 


0. 


06 


25.8 


5.7 


3. 


02485 





.00006 


1 


W2-0Q 6 


3888.38 


0. 


03 


83.9 


4.4 


3. 


02487 





.00003 




W3-0Q 6 


3813.22 


0. 


05 


55.7 


7.2 


3. 


02484 


0. 


.00005 




Ll-OP 6 


4403.42 


{) 


04 


25.4 


4.1 


3 


02486 


() 


00004 




L3-0R b 


4280.27 


0. 


.03 


65.6 


3.8 


3 


02485 


o 


00002 




L3-0P a,b 


4984 8Q 


0. 


04 


24.2 


4.2 


3 


02486 


o 


00004 




L7-0R b 


4078.95 


0. 


03 


76.4 


4.8 


3 


.02485 


0. 


.00003 




L10-0P b 


3955.77 





04 


46.8 


5.0 


3. 


02484 





.00004 




L14-0P a ' b 


3813.61 


0. 


09 


24.6 


7.2 


3. 


02484 





.00009 


2 


W0-0R a ' b 


4061.19 


0. 


.05 


34.4 


5.5 


3. 


02488 





.00005 




W0-0Q a ' b 


4068.87 





04 


50.7 


6.7 


3. 


.02484 





.00004 




W0-0P 6 


4073.85 


0. 


09 


10.3 


4.1 


3. 


.02488 


0. 


.00009 




Wl-OQ 6 


3976.45 





04 


36.6 


4.4 


3. 


02487 


0. 


.00004 




L3-0R a ' b 


4286.44 


0. 


06 


14.9 


4.7 


3. 


.02485 


0. 


.00006 




L3-0P 6 


4294.12 


0. 


07 


12.6 


3.9 


3. 


.02485 





.00007 




L4-0R a ' b 


4232.14 


0. 


.05 


24.9 


5.0 


3. 


02487 


0. 


.00005 




L5-0R a ' b 


4180.57 


{) 


04 


32.2 


4.4 


3. 


.02485 





.00004 


3 


W0-0Q 6 


4075.88 


{) 


03 


65.7 


5.0 


3. 


02485 





.00003 




L2-0R b 


4353.77 


0. 


03 


36.1 


4.6 


3. 


02488 


0. 


.00003 




L2-0P a 


4365.17 


0. 


07 


8.6 


5.7 


3. 


02481 





.00007 




L3-OR b 


4296.43 


0. 


03 


51.4 


4.6 


3. 


.02484 


0. 


.00003 




L3-0P a ' b 


4307.18 


0. 


05 


22.0 


4.5 


3 


02487 


0. 


.00005 




L4-0R b 


4252.15 


0. 


03 


61.1 


5.4 


3. 


02486 





.00003 




L4-0P 


4252.14 





09 


8.4 


6.0 


3. 


02486 





.00008 




L6-0R b 


4141.54 


0. 


03 


60.3 


4.6 


3. 


.02487 





.00003 




L7-0P 6 


4103.35 


0. 


05 


32.2 


4.0 


3. 


02484 





.00005 




L8-0P a ' b 


4058.61 


0. 


05 


30.1 


5.5 


3 


02484 





.00005 




L12-0R a ' b 


3894.77 


0. 


04 


49.0 


6.2 


3. 


02488 





.00005 




L15-OR a ' b 


3795.31 


0. 


08 


20.8 


5.0 


3. 


02486 


0. 


.00008 


4 


L4-0P a ' b 


4268.64 


{) 


05 


24.1 


5.9 


3. 


02482 





.00005 




L8-0P a 


4074.21 


0. 


09 


8.6 


5.0 


3. 


02484 


0. 


.00009 


5 


Wl-OQ a ' b 


4004.44 


0. 


.06 


25.7 


5.4 


3. 


.02489 


0. 


.00006 




L7-0P a ' b 


4138.55 





.10 


8.5 


4.7 


3. 


02484 


0. 


.00010 




L12-0R a ' b 


3923.76 


{) 


08 


13.9 


4.4 


3. 


02484 





.00008 



Note. — a Lines are used in the Voigt-fitting procedure to estimate the broadening b- 
parameter. b Lines are used in Figure 3. 
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Table 2 

Molecular hydrogen and metal abundance measurements in DLAs 



# 


QSO 


z. 


lbs 


[Cr/Zn] 


[Fo/Zn] 


log Af(Hi) 


logJV(H 2 ) 


log /h 2 


1 


Q1328+307 





.692 


-0.44 ± 0.10 a 




21.28 ± 0.10 a 


< 15.70 b 


< -5.28 


2 


Q0454+039 





.860 


0.14 ± 0.12 c 


0.01 ± 0.12 c 


20.69 ± 0.06 c 


< 16.38 b 


< -4.01 


3 


Q0935+417 


1 


.373 


-0.10 ± 0.16 d 


-0.27 ± 0.15 d 


20.40 ± 0.09 d 


< 16.40 b 


< -3.70 


4 


Q1157+014 


1 


.944 




-0.42 ± 0.11 s 


21.80 ± O.lO-f 


< 14.71" 


< -6.79 


5 


Q0013-004 


1 


.973 


< -0.95 s 


-1.10 ± 0.13 s ' h 


20.70 ± 0.05* 


19.84 ± 0.10* 


-0.66 ± 0.10 


6 


Q0458-020 


2 


.040 


-0.36 ± 0.05 j 


-0.46 ± 0.07 j 


21.65 ± 0.09 s 


< 15.90 b 


< -5.45 


7 


Q0100+130 


2 


.309 


-0.12 ± 0.05 j 


-0.22 ± 0.06 j 


21.40 ± 0.05 fc 


< 15.70 b 


< -5.40 


8 


Q1232+082 


2 


.338 




-0.86 ± 0.13'' m 


20.90 ± 0.10" 1 


17.18 ± 0.10" 1 


-3.42 ± 0.14 


9 


Q0841 + 129 


2 


.374 


-0.12 ± 0.08 e ' j 


-0.29 ± 0.08° 


20.95 ± 0.09 j 


< 14.67" 


< -5.98 


10 


Q0112+029 


2 


.423 


-0.75 ± 0.18 s 




20.95 ± 0.10 s 


< 16.00 b 


< -4.65 


11 


Q1223+178 


2 


.466 


-0.27 ± 0.14 s ' e 


-0.21 ± 0.11" 


21.52 ± 0.10 e 


< 14.30" 


< -6.92 


12 


Q0528-250 


2 


,811 


-0.46 ± 0.14" 


-0.46 ± 0.12"'° 


21.35 ± 0.10"' p 


16.77 ± 0.09°'° 


-4.28 ± 0.13 


13 


Q0347-382 


3 


.025 


-0.37 ± 0.06" 


-0.19 ± 0.06" 


20.626 ± 0.005" 


14.613 ± 0.022" 


-5.71 ± 0.02 


14 


Q0000-262 


3 


,390 


0.06 ± 0.07* 


0.03 ± 0.07* 


21.41 ± 0.08^ 


13.94 ± 0.06"'™ 


-7.17 ± 0.13 



Note. — Column densities N(Hl) and 7V(H2) arc given in cm 2 . 
References: a Meyer fc York (1992); b Gc & Bcchtold (1999); c Pcttini et al. (2000); d Meyer et al. (1995); e Pctitjean et al. 
(2000); / Wolfe & Briggs (1981); 9 Pettini et al. (1994); h Centurion et al. (2000); * Gc & Bcchtold (1997); j Prochaska & 
Wolfe (1999); k Wolfe et al. (1994); ' Ge et al. (2001); m Srianand et al. (2001); " Lu et al. (1996); ° Srianand & Petitjean 
(1998); v M0ller & Warren (1993); s Levshakov & Varshalovich (1985); " this paper; * Molaro et al. (2000); " Levshakov 
& Foltz (1988); " Levshakov et al. (2000); *" Levshakov et al. (2001) 
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Table 3 

Column densities and relative abundances in the H 2 -bearing cloud at z abs = 3.025 toward 

Q0347-3819 



Ion 


log N, cm 2 


(X/H) 6 


[X/H] 


[X/Zn] 


Hi 


20.402 ± 0.006° 








Ci 


< 11. 6 t4 








Cn* 


13.354 ± 0.023 








Cn 


15.703 ± 0.024 


-3.48 ±0.06 


-1.22 ± 0.06 


0.29 ± 0.08 


Ni 


14.664 ± 0.007 


-4.08 ± 0.06 


-1.66 ± 0.06 


-0.15 ± 0.08 


Oi* 


< 11.9 1 '* 








Oi 


16.410 ± 0.010 


-3.17 ± 0.06 


-0.82 ± 0.06 


0.69 ± 0.08 


Sin* 


< 12. 4 f '* 








Sin 


15.015 ± 0.013 


-4.44 ± 0.05 


-0.95 ± 0.02 


0.56 ± 0.06 


Pn 


12.779 ± 0.032 


-6.44 ± 0.06 


-1.18 ± 0.07 


0.32 ± 0.09 


Ar I 


13.760 ± 0.020 


-5.48 ± 0.04 c 


-1.16 ± 0.05 


0.35 ± 0.07 


Cm 


11.698 ± 0.046 


-6.31 ± 0.01 


-2.39 ± 0.05 


-0.89 ± 0.07 


Foil 


13.690 ±0.032 


-4.50 ± 0.01 


-2.21 ± 0.03 


-0.71 ± 0.06 


Znn 


11.565 ± 0.037 


-7.33 ± 0.04 


-1.51 ± 0.05 





Note. — "Indicated are the internal errors, t la level; * based on the Keck/HIRES data obtained by Prochaska & Wolfe (1999). 
References : b Grevessc & Sauval (1998); c from Sofia fc Jenkins (1998). 
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Table 4 

Total column densities and relative abundances at z abs = 3.025 toward Q0347-3819 t 



Ion 


log N, cm 2 


(X/H)& 


[X/H] 


[X/Zn] 


[X/H] 6 


[x/ar 


[X/H] = 


Hi 


20.626 ± 0.005 x 














Cii 


15.929 ± 0.024 


-3.48 ± 0.06 


-1.22 ± 0.06 


0.29 ± 0.08 




> -2.4 


> -2.2 


Cn* 


13.581 ± 0.024 














Ni 


14.890 ± 0.007 


-4.08 ± 0.06 


-1.66 ± 0.06 


-0.15 ± 0.08 


-2.1 


-2.5 




Oi 


16.636 ± 0.011 


-3.17 ± 0.06 


-0.82 ± 0.06 


0.68 ± 0.08 


> -1.5 


> -2.4 


> -2.2 


Sin 


15.240 ± 0.013 


-4.44 ± 0.05 


-0.95 ± 0.02 


0.56 ± 0.06 


-1.2 


> -1.9 


-1.4 


Pn 


13.005 ± 0.031 


-6.44 ± 0.06 


-1.18 ± 0.07 


0.32 ± 0.09 








An 


13.986 ± 0.020 


-5.48 ± 0.04 a 


-1.16 ± 0.04 


0.34 ± 0.07 








Cm 


12.441 ± 0.035 


-6.31 ± 0.01 


-1.88 ± 0.04 


-0.37 ± 0.06 








Fen 


14.432 ± 0.009 


-4.50 ± 0.01 


-1.69 ± 0.01 


-0.19 ± 0.06 


-1.9 


> -2.0 


-1.7 


Znn 


11.792 ± 0.036 


-7.33 ± 0.04 


-1.50 ± 0.05 




< -0.8 


< -1.2 





Note. — * Listed values were measured in the range —50 km s 1 < v < 20 km s 1 (sec Figure 10). x Indicated 
are the internal errors. References : a from Sofia & Jenkins (1998); b Centurion et al. (1998); c Ledoux ct al. 
(1998); d Grcvcssc & Sauval (1998); e Prochaska &c Wolfe (1999) 
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Table 5 

Transition probabilities and excitation rate coefficients for AJ = 2 transitions 

BETWEEN THE FIRST FIVE ROTATIONAL STATES 



J 




„T=400K 1 cm 3 "\ b 
9j^J+2 \~s~ ) 


T=1000K / cm 3 "\ b 

9j^j+2 ^— —) 


T=1000K / cm 3 "\ c 
9.T^.T + 2 \~s~ ) 


T = 1000K ( cm 3 "\ d 
</.T^.T + 2 \~s~ ) 





2.94(-ll) 


1.14(-11) 


6.44(-ll) 


4.82(-ll) 


1.27(-10) 


1 


4.76(-10) 


1.79(-12) 


1.60(-11) 




4.07(-ll) 


2 


2. 76 (-9) 


5.16(-13) 


7.35(-12) 


6.70(-12) 


1.78(-11) 


3 


9.84(-9) 


1.76(-13) 


4.02(-12) 




8.33(-12) 



Note. — The negative numbers in parenthesis arc the powets of ten. References : a Turner ct al. (1977); this 
paper; c Elitzur fc Watson (1978); d Nishimura (1968) 
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Table 6 

Comparison of the mean values and the standard deviations estimated by the Monte 
Carlo (MC) and Ay 2 methods in case of correlated measurements 



parameter, x 




(i)mc 


CMC 




a ±x 2 


6*, km s- 1 




3.04 


0.29 


2.80 


0.45 


JV(J=0), cm" 


2 


1.90{ 3 


0.29i3 


1.90i3 


0.35i3 


JV(J=1), cm" 


2 


1.84 14 


0.13i4 


I.6O14 


O.I814 


JV(J=2), cm" 


2 


5.39i 3 


0.28i 3 


5.40i3 


0.35i 3 


JV(J=3), cm' 


2 


1.24 14 


0.07 14 


1.22i4 


0.05i4 


iV(J=4), cm" 


2 


2.35i3 


0.44i3 


2.6O13 


0.43i3 


JV(J=5), cm" 


2 


2.90i3 


0.53i3 


2.90i3 


0.70i3 



Note. — * The H2 column densities at different J levels and 
the corresponding broadening b parameter are listed in col- 
umn (1). t Values like ' at mean a X 10 . 
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Wavelength (A) 
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0.96 
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• - 
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X 
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• 

la 
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1 
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0.92 
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. • 1 
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i 1 i i 

•■ 


X 


• 

• 
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